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Abstract The age of an individual star cannot be measured, only estimated through mostly 
model-dependent or empirical methods, and no single method works well for a broad range of 
stellar types or for a full range in age. This review presents a summary of the available techniques 
for age-dating stars and ensembles of stars, their realms of applicability, and their strengths and 
weaknesses. My emphasis is on low-mass stars because they are present from all epochs of star 
formation in the Galaxy and because they present both special opportunities and problems. The 
ages of open clusters are important for understanding the limitations of stellar models and for 
calibrating empirical age indicators. For individual stars, a hierarchy of quality for the available 
age-dating methods is described. Although our present ability to determine the ages of even 
the nearest stars is mediocre, the next few years hold great promise as asteroseismology probes 
beyond stellar surfaces and starts to provide precise interior properties of stars and as models 
continue to improve when stressed by better observations. 
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1 WHY AGES MATTER 

We can directly measure many key physical properties of a star. The masses of 
stars in binary orbits can be determined from basic physics, and we can infer the 
mass of most single stars to ~ 10% just from a spectral type, particularly for 
the main sequence. Determining the composition of a star is not easy, but it is 
straightforward, and the process and its limitations are well understood. But age, 
the third key determinant of a star's physical state, is another matter all together. 
The Vogt-Russell theorem tells us that the physical state of a star results from its 
mass and composition. Other factors (rotation, companionship, magnetic fields) 
matter too, but mass and composition dominate. The composition of a star 
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changes with time, and so age influences the state of a star, but less directly. 
Time is a medium in which the star inexorably changes, but age is not the direct 
agent of that change. 

All of this means that a star's age cannot be measured, it can only be estimated 
or inferred. This is underscored by the fact that we have exactly one stellar age 
that is both precise and accurate, and that, of course, is for the Sun. But the 
Sun itself does not reveal its age; it is only because we can study solar system 
material in the laboratory that we can limit the Sun's age. We can do that for 

no other star. Sidebar 1: Why Care 

Yet age lies at the heart of stellar evolution, for how can there be evolution About Stellar Ages? 
without the passage of time? Sometimes our desire to know ages is satisfied if we (see p. I46p 
can place events in the correct sequence and understand the relative duration of 
phases. But often we need to tie together observations of independent objects to 
arrive at an understanding of a process. Some examples include: 

• The formation and evolution of proto-planetary disks appear to occur in 
the first ~ 100 Myr of a star's life, with debris disks forming later. At 
present we can just barely limit this time-scale with the methods available, 
but clearly one would like to be able to see differences in the formation 
processes and time-scales from place to place and star to star if we are to 
understand fully the physics of the processes involved. 

• Many gas-giant planets have now been found around nearby solar-type 
stars. We can estimate rough ages for those stars, but a better under- 
standing of, for instance, the dynamics of such systems and the inward 
migration of planets requires accurate ages. 

• We could better understand the behavior and future of our Sun if we could 
create a true cohort based on mass, composition, and age. At present, 
comparisons of the Sun to other stars often rely on the very stellar properties 
(such as activity, rotation, or lithium abundance) that one is trying to 
understand; the process is self-referential. 

• To fully understand the Galaxy's history of formation, enrichment, and 
dynamics, we need to be able to assign ages to individual stars in its com- 
ponents. Low-mass stars live long enough to be present from all epochs of 
star formation, making them attractive for these purposes. Recent studies 
have applied some of the techniques discussed in this review to Galactic 
problems, with differing and sometimes conflicting results. 

• The search for life and any understanding gained from observing signs of 
life beyond our solar system has a profound need to measure stellar ages if 
we hope to gain insights into biological evolution. 

These examples come to mind because of my own interests, and the reader 
can supply many more. A stellar chronologist will not lack for work. A full 
understanding of many phenomena in astrophysics means at least getting the 
sequence of events right, and, preferably, knowing time-scales. Many problems 
are connected to individual stars or the nearest stars, and we need reliable ways to 
get at the ages of such objects. I also note that the overall cosmic age scale, which 
is most unambiguously constrained by the age of the solar system, is perhaps the 
single most contentious point in all of science for those who are science deniers, 
and an understanding of that datum and of ages in general is important for us 
as involved citizens. 
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This review will discuss a broad range of stars, both as individuals, and as 
clusters, associations, and groups. Cluster ages are important in themselves, of 
course, but here I use clusters mostly as benchmarks to calibrate other age estima- 
tors. I will emphasize late-type stars because they are numerous, astrophysically 
interesting in themselves, and are present from all epochs of star formation in the 

Galaxy. Sidebar 2: A Challenge 

for Stellar Astrophysics 
1.1 What Do We Mean By Age? (see p. ST]) 

When considered in detail, the age of a star is inherently uncertain if only because 
the point at which its age begins is ill-defined. Theory might use the point when 
hydrostatic equilibrium is established, but that is not easy to pinpoint observa- 
tionally. One baseline that is reasonably well-defined in th e observations of pre- 
main sequence stars is the deuterium-burning "birthline" of lStahler at the 



point where proto-stars lose the ir shrouds and reveal themselve s. At that point a 
1.0 Mq star is ^ 10^ years old. Wuchterl &: Tscharnuter (20031 ) have proposed a 



time zero-point for when an object's photospher e first exists, which they take to 
be when the Rosseland mean opacity is 2/3. The lWuchterl fc Tscharnuter f 20031 ) 



zero-point precedes the birthline by perhaps ^ 10^ years and is dependent on the 
star's accretion history ( Froebrich et al. 20061 ). These differences are critical for 



fully understanding the earliest phases of star and planet formation, but that de- 
gree of uncertainty becomes essentially meaningless for the stars to be discussed 
here, which are ~ 1 Myr old and more. 

Some have suggested using the arrival of a star on the zero-age main sequence 
to define = 0." That is a well-defined point in a star's life, both theoretically 
and observationally, but it is midway through the star's life and would require 
the use of negative ages. We cannot yet define well the true birth of a star in 
a dense cloud, but the birthline is reasonably well-defined, both in models and 
observationally. Starting with t = at that point makes it convenient to work on 
a logarithmic scale. As we will see, most of the various age estimators are also 
best used logarithmically. 

The concept of "age" is different for white dwarfs (WDs) because the age is 
based on the cooling of the WD (see Sec. I9.2.2p . This calculation is based 
on well-understood physics, which is what makes WD cooling ages attractive, 
but the cooling does not start until the WD forms, and the time elapsed up to 
that point depends on more conventional stellar modeling that is subject to a 
number of uncertainties. This is not important for older clusters, where the WD 
progenitors are massive stars for which the pre- WD phase is a small fraction of 
the cluster's age. 

1.2 Goals of This Review 

The need for reliable stellar ages has always existed, but recent work has added 
significantly to the observations available to address these problems, and the 
quality of the models has kept pace. The problems are not solved, but we probably 
have a better understanding of the limitations. Also, as I note at the end, I expect 
major advances in this area over the next few years. 

I can not review the ful l history of this subject, but a key motivating paper was 
that of Skumanich f 19721 ) ■ who looked at some age-related observable properties 



of solar-type stars. He had little data to work with at the time, but the power-law 
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relations he plotted for rotation and activity (and an exponential for depletion of 
lithium) instigated much subsequent work. 

Despite the unambiguous importance of knowing the ages of stars, the subject 
has been the focus of only a very few topical symposia and few or no previous 
general reviews. I will start with an outline of the available techniques that can be 
applied to individual stars and a scheme for understanding them. I will emphasize 
"Population I" or "thin disk" stars and low- to intermediate mass stars simply 
because they predominate in the solar neighborhood and are central to many of 
the questions for which an age is part of the answer. I hope to improve on the 
status quo by compiling the available information in one place and by attempting 
to systematically lay out alternatives and their advantages and shortcomings. I 
will compare the different methods when possible, and I will point out areas that 
can benefit from attention. 

This review will focus mostly on late-type dwarfs and pre-main sequence stars. 
Elderly individual stars will be discussed, but the roader topic of Population II 
stars is beyond the scope of this. Similarly, white dwarfs (WDs) are also very use- 
ful for age-dating clusters and populations us ing fundament al ly different physic s, 
but they have been well reviewed elsew here ( Koester (20021 ) . Ivon Hippel (20051 ) ■ 
Winget Kepler r2008l l. ISalaris r2009l l) and so will only be mentioned here (see 



also Sec. I9.2.2p . Evolved stars present some special problems of their own and 
also will not be discussed here. However, more on the ages of these objects (ari d 
on stellar ages in general) can be found in Mamaiek. Soderblom &: Wvse (20091 ). 
which contains review papers from lAU Symposium 258, "The Ages of Stars." 

Evolved stars, because of their high luminosities, are often used as tracers of popu- 

lations, and a discussion of the ages of populations can be found in lTolstoy. Hill &: Tosi (20091 ). I 
Non-stellar objects will not be treated here. 

I will also briefly re-examine open clusters and cluster ages as the benchmarks of 
other methods. The correct age scale to be used for OCs remains fundamentally 
uncertain to a significant degree, as we will see, but for the most part other 
age estimators are at least satisfactory if they provide results on a scale that is 
consistent with the OCs. Other ways to calibrate age indicators are also discussed. 



1.3 Some Helpful Terms and Abbreviations 

A few abbreviations are used here so often that this is a good place to define 
them. ZAMS is, of course, the zero-age main sequence, the point in a star's life 
when hydrogen starts to fuse into helium in the core. Pre-main sequence (PMS) 
stars are young objects that have not yet reached the ZAMS, and main sequence 
(MS) stars are more generically ZAMS and later. We distinguish between the 
Hertzsprung- Russell diagram (HRD), which is in physical units of luminosity and 
temperature, and the color-magnitude diagram (CMD), which is in observational 
units. One is better suited to some discussions than the other, and the transfor- 
mations between the two are not trivial. A frequently-discussed locus in either 
the HRD or CMD is the main-sequence turnoff (MSTO), the region in an OCs 
HRD where stars stop getting hotter and bluer and an isochrone bends back to 
the red for more massive stars. We will use r and log r to denote an age, in years. 

Open clusters (OCs) are discussed because they are so important as age bench- 
marks. OCs are well-defined groups of stars that are clearly in close proximity 
to one another and which are defined observationally by having only a small 
dispersion in radial velocity and proper motion. OCs may not be physically 
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bound objects, but they form an easily-recognized ensemble. Less well-defined 
are groups and associations, terms which are used interchangeably here. The 
youngest stars we know of are found in star-forming regions, in close physical 
proximity to one another, but those groups are clearly too sparse to be bound. 
At somewhat greater ages (~ 10 Myr), loose an d very sparse groupings h ave re- 
cently been identified in the solar neighborhood (jZuckerman &: Song 20041 ). Most 



of the bona fide OCs have ages from 10 to 200 Myr, older clusters being 
rare due to their tidal disruption in the Galaxy. These sparse groupings are at 
least as much the subject of study for their astrophysical properties as they are 
benchmarks for age indicators. 

2 GENERAL COMMENTS 

2.1 Attributes of a Useful Age Indicator 

We cannot n ieasure stellar ages directly and so we use other tactics to address 
the problem. Barnes (20071 ) has discussed some of the attributes of a measurable 



quantity that would make it an ideal age indicator. They include having a well- 
defined quantity that can be measured for single stars and that is sensitive to 
age but not other factors. The quantity should be able to be calibrated and 
with known uncertainties. The relation between the indicator and age should be 
invertible so that age can be calculated, and when applied to a coeval sample 
one should get consistent ages. These criteria are all valid but address an ideal, 
and to that list should be added that the physics of the relation between the 
indicator and age is understood and that establishing the relation involves the 
fewest possible assumptions. 

There is no single ideal age indicator, and for real stars in the real Universe 
we have to settle for less, often much less. Our errors may be dominated by 
systematic effects since we do not usually know what inherent physical relation 
to expect. The sensitivity of the indicator may be useful for only a limited age 
range or mass range. Inherent scatter in the quantity measured may mean that 
the age indicator is not very useful for single stars but can work if at least a 
few stars are being considered as an ensemble. Clusters are ordinarily used to 
calibrate and to test how well a coeval population behaves, but open clusters 
more than ~ 1 Gyr old are inherently rare and tend to be distant. 

We may be willing to get by with much less. In many cases it suffices to get 
the ordering of ages right, even if the age scale remains uncertain. In other cases 
it may be enough to set limits to an age. There are no perfect age indicators, 
and each decade of stellar age and each stellar type presents its own problems 
and has its own relevant indicators. In order of usefulness, age indicators can 
be used, first to classify or associate, as is done, for instance, when a strong Li 
feature is used to confirm that a star is a member of the T Tauri class. Second, 
an indicator can be used to limit an age, perhaps by imposing an upper bound. 
Next, it can be helpful to sort ages in order to get the order right. Finally, we 
can hope to quantify ages so as to specify a numeric age or age range. 

2.2 A Scheme for the Quahty of Age Estimation Methods 

I break down the available age estimation methods into five levels of quality, 
based on how straightforward it is to go from the measurement to the age and 
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how well we understand the physics at each step. Each of the methods listed will 
be explicated further. The table below summarizes my own view of how suitable 
the different methods are within different domains as a means of illustrating how 
the subjects will be presented. 

2.2.1 FUNDAMENTAL AND SEMI-FUNDAMENTAL AGES An age 
is fundamental if the underlying physical processes are completely understood and 
all the needed observations are available. The one and only fundamental age we 
have to work with is that of the Sun (see sidebar) , and it is the result of measuring 
the decay products of long-lived isotopes in solar system material, something we 
can do for no other star. 

There are two stellar age estimation methods that can be considered semi- 
fundamental because they involve making only a few assumptions, those assump- 
tions appear to be well-founded, and the assumptions do not influence the derived 
ages to a large degree. The first, nucleocosmochronometry, applies to individual 
stars of the thick disk and halo and involves measurement of the decay of U or, 
especially, Th. This is done for individual stars. The method has been used only 
on metal-poor stars because it is only in those cases that the weak Th and U 
absorption lines can be measured in the presence of many blends. The complex 
spectra of stars in and around the Th and U lines remains an important limitation 
on this method. 

The second semi-fundamental technique can only be applied to a group of 
young stars and involves tracing the Galactic orbits of the stars back in time. An 
assumption must be made about the Galactic potential through which the stars 
move, but that is not critical. The age of the group is assumed to be when the 
stars were in closest physical proximity to one another. This method of expansion 
ages works only for young groups because massive objects in the Galactic disk 
tidally disrupt groups and clusters on a time-scale of ~ 200 Myr, leading to "disk 
heating." If a group has an age less than ~ 50 Myr we can be confident it has not 
yet been so disrupted. Also, the uncertainties involved in calculating the stellar 
kinematics make it difficult to go back farther in time than 30 Myr. 

2.2.2 MODEL-DEPENDENT AGES The majority of the age estimation 
techniques used have built-in model dependencies. The best-known method for 
clusters - isochrone fitting - starts with a series of calculated models and varies 
the properties of those models to achieve the best overall match to the observa- 
tions. In addition to the quality of the models themselves, there is also generally 
a calibration dependency in that T^s, for instance, must be related to an observed 
color or colors, and reddening taken into account since interstellar extinction is 
not a property of the models themselves. Every step of the process adds a link 
to the chain. Also, the process of isochrone fitting is itself a test of the mod- 
els being used and our understanding of the underlying physics. Examples of 
model-dependent methods include: 

• Isochrone fitting, for clusters, or isochrone placement, for individual stars. 
As we will describe below, for clusters isochrone fitting has the advantage 
of using a cluster's entire HRD, particularly key features that depend on 
specific physical processes. This is what makes the HRD such a valuable 
tool for understanding stellar evolution, of course, and for testing mod- 
els of stellar physics. However, the quality of the isochrones depends on 
the completeness and accuracy of the models, and missing or inadequately 
understood processes can lead to systematic error. 



Table of methods: see 
separate PDF 
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• Asteroseismology - the detection of oscillation modes in stars - is espe- 
cially promising for single solar-type stars and for older stars because the 
low-order modes one can detect in an unresolved object pass through or 
near the star's center, making them an indicator of the star's central den- 
sity, which is to say its age. Interpreting asteroseismological observations 
requires detailed models, but the physics in those models is generally differ- 
ent than for MSTO stars in clusters and is well-understood, being for stars 
similar to the Sun. 

• The lithium depletion boundary (LDB) seen in very young open clusters is 
the point at the bottom of the main sequence where Li reappears. Stars 
more massive than the LDB astrate all their Li in approaching the MS, while 
less-massive objects cannot reach the necessary internal temperature. The 
observations needed to locate the LDB in a CMD are challenging because 
such objects are inherently very faint, but the underlying physics is fairly 
simple and well-understood. 

2.2.3 EMPIRICAL AGES Empirical ages are based on observed relations 
between a measured property and age. The full physics of the relation is not 
understood and so there is no predictive power, but there are at least reasonable 
and plausible scenarios to explain what is seen. Empirical ages must be empiri- 
cally calibrated by measuring the relevant property in members of open clusters. 
Because of the scarcity of old open clusters, combined with the difficulty of mea- 
suring the various properties in older stars, the Sun is often used to establish the 

relation for older stars and this adds an additional ass umption. 

The most direct of these empirical methods is what iBarnes has called 



gyrochronology: using rotation periods to derive an age. Less direct are methods 
using various forms of stellar activity (seen in the cores of the Ca ii H and K 
or the Mg II h and k lines, in Ha, and in x-rays, primarily). Using activity 
is less direct because the activity appears to be dependent on the rotation. In 
addition, activity diagnostics are inherently variable. The last empirical method 
uses lithium abundances in F, G, and K stars. 
2.2.4 STATISTICAL AGES Two statistical correlations have been dis- 
cussed extensively: the age-metallicity relation (AMR) and the net increase in 
Galactic space motions with time ("disk heating"). It is not clear than an AMR 
relation actually exists in that the observed correlation may be due to an over- 
all connection between Galactocentric radius and metallicity, with stars getting 
scattered into the solar neighborhood from the many places in our Galaxy where 
they formed. At any rate, the correlation is at best a coarse one and is only 
suggestive for an individual star. As a counterexample, there are recently-formed 
stars with metallicities about a factor of two below solar. 

There is also a general trend of stellar space motions increasing with age, but 
note that both the Sun and the a Centauri system are ~ 4 — 5 Gyr old yet have 
space motions near the Local Standard of Rest. Again, a high space motion is 
suggestive of an old age, but it is also possible for individual younger stars to 
have been ejected from a cluster and for older stars to appear unaffected by these 
forces. However, the net kinematic properties of, for instance, a volume-limited 
sample or another well-defined ensemble can be used to estimate a mean age. 
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3 FUNDAMENTAL AND SEMI-FUNDAMENTAL AGES 

As noted, the one and only fundamental age in stellar astrophysics is that of 
the Sun, based on measuring abundances and daughter products of radioactive 
nuclides in meteoritic material. Similar nuclide abundances can be measured 
in the spectra of stars, but without full knowledge of all the nuclides involved 
some critical assumptions must be made. Those assumptions appear to be rea- 
sonably sound on both theoretical and practical grounds, leading to the first 
semi-fundamental technique, nucleocosmochronometry. 



3.1 Nucleocosmochronometry 

The ability to detect and measure Th and U in metal-poor stars has made nu- 
cleocosmochronometry particularly attractive for studying stars formed in the 
earliest epochs of the Galaxy and the processes that have enriched them. Such 
stars are generally much too far away to have measured parallaxes that could 
lead to luminosity determinations and isochrone placement, but are compelling 
for understanding the time-scales of nucleosynthesis in the early Galaxy and for 
comparing what we believe to be the oldest stars to the ages of globular clusters 
(|VandenBerg. Bolte Stetson 199^ ). 



Nucleocosmochronometry derives stellar ages by measuring the decay of long- 
lived isotopes. One seeks an isotope with a half-life comparable to the age of 
the object under study so that a measurable quantity remains; the isotope must 
offer detectable features; and the element's abundance should be dominated by 
that isotope. The decay of these isotopes involves fully understood physical pro- 
cesses. However, an initial abundance for the elements must be assumed because 
it cannot be measured directly. This is ordinarily done by scaling from other 
r-process abundances and assuming production ratios of the elements being mea- 
sured. There remains significant uncertainty in these steps. Different calculations 
c an lead to production ratios that diffe r by nearly a factor of two. 



Cowan. Thielemann &: Truran (19911 ) provided a thorough review of th e his 



tory and early applic ation of this method, and Beers &: Christlieb f2005l l and 



Frebel &: Kratz (20091 ) provide additional information. Ludwig et al. (20091 ) present| 



a summary of recent work on nucleocosmochronometry and the errors associated 
with it. The two species used are ^^^U (ri/2 = 4.47 Gyr) and '^^'^Th (ri/2 = 14.05 
Gyr). Of the two, Th is more easily observed, with a Th II feature at 4019 A. 
U is more d ifficult to detect e xcept in some cases of highly-enhanced r-process 
abundances ( Frebel et al. 20071 ). Other nuclides with suitable half-lives have been 



sought that present good spectrum lines but none have been found. 

The U and Th isotope abundances are determined relative to other rare earths 
such as Eu, Os, Hf, or Ir. Ideally the comparison element is one formed in 
just the same way as the T h, say, but this con nection is an assumption that 
adds uncertainty. Studies by Cowan et al. (19991 ) and others cited there tend to 



support the idea that r-process elements are formed in consistent proportions. 
Kratz et al. !(2003) argue th at Hf should work b est because it is so close in mass 



to Th, yet the analysis of iLudwig et al. does not bear that out, with 



Th/Hf ratio s producing consiste ntly unphysical results. An illustrative case is 
provided by ISneden et al. !(2003) in their reanalysis of the very metal-poor star 



CS 22892-052. They derived a net age of 12.8 ± 3 Gyr. Thi s is, of course, consis - 
tent with WMAP's age for the Universe of 13.7 ± 0.2 Gyr (jBennett et al. 20031 ). 
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but the individual ages from various ratios in ISneden et al. (200i ) differed by 
nearly a factor of two (from 10.5 to 19.2 Gyr). Another very metal-poor star 
with enhance d r-process elemeii ts, HE 1523-0901, presents a better case, having 
detectable U (iFrebel et al. 200?! ) and highly consistent ages from both U and Th, 



with an average of 13.2 Gyr. 

The other examples of recent analyses cited in iLudwig et al. present 



much less consistency. Often the calculated ages exceed 20 Gyr and only in a 
few cases are they below the WMAP age. In a few cases the derived ages are 
~ 1 Gyr or less or even formally negative. No one comparison element by itself 
provides better consistency. The promise of nucleocosmochronometry to provide 
ages for, say, globular clusters that are independent of isochrone fitting will be 
fulfilled only with significant additional effort. Systematic errors from different 
calculated production ratios should be able to be reduced, and the improving art 
of spectrum synth esis will lead to better results as well. For comparison, about 
a half century ago Fowler &: Hovle (1960l ) determined an age for the Galaxy of 



15 lb 3 Gyr from nuclear analysis of solar system material. 

In principle nucleocosmochronometry would be ideal for old main sequence 
stars, where few other age estimation techniques are useful, but in practice this 
method can generally be applied only to metal-poor stars because of the difficulty 
of measuring the U or Th feature in solar-metallicity stars in the presence of the 
many blending features. In addition, the r-process elements in a star like the 
Sun will have come from many separate events, while in the oldest stars formed 
early in the Galaxy's history fewer nucleosynthesis events will have taken place 
before the star formed. Also, in some cases the comparison element will have a 
contribution from the s-process as well as the r-pro cess. Nevertheless, t he Sun's 
age has been calc ulated as a test of t he procedures. iLudwig et al. ilooi) use the 



measurements of ICaffau et al. (20081 ) to derive solar ages from 1.7 to 22.3 Gyr. 



The Th/Eu ratio appears to yield the most reasonable results, but with poor 
precision. 

Th/Eu ratios have also been used by del Peloso. da Silva &: Aranv-Prado f2005l ).| 



del Peloso. da Silva fc Porto de Mello (200,51 ). and ldel Peloso et al. f 200,51 ^ to de- 



termine an age for the Galactic disk of 9 it 2 Gyr. In this case the age is not 
from individual stars but from matching the observations to models of Galactic 
chemical evolution that take several effects into account that enrich and destroy 
Th and Eu with time. 

3.1.1 SUMMARY OF NUCLEOCOSMOCHRONOMETRY 

-|- The method is applicable to the Galaxy's oldest stars, independent of distance. 
-|- The method is applicable to individual stars. 
-|- The physics of the process is well-understood. 

— Spectra of high resolution and high signal-to-noise are required to detect the 

weak features, limiting the method to brighter stars. Generally only a single 
feature of an element is detected. 

— The method works only for low-metallicity stars because of line blending. Very 

few stars have been found with detectable Th or U features. 

— The uncertainties are at least 20%. iLudwig et al. r2009l ) estimate that with 



even favorable errors that the uncertainty in age for an individual star is 
about 2.5 Gyr. 
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— Significant systematic errors remain, particularly in the expected production 
ratios of rare earth elements. 

3.2 Kinematic or Expansion Ages 

Kinematic ages are determined by tracing the motions of a group of stars into the 
past to determine when they were in closest physical proximity, which is assumed 
to be the time of formation. Alternatively, one can plot proper motion in decli- 
nation, say, against declination to get a linear expansion rate, which then implies 
an age. The concept goes back to work on the exp ansion of OB associations 
by Ambartsumian and Blaauw around 1950 (see, e.g., Blaauw (19641 )) . although 



the early results were inconsistent and sometimes in conflict with evolutionary 
ages for massive stars. A related age-dating method connects runaway stars back 
to their point of formation and assumes linear motion. That topic will not be 
discussed here as runaway stars are important in their own right as tracers of 
dynamical processes. 

The attraction of the kinematic method, which obviously can only be applied 
to a group of stars, is that it is independent of any stellar physics since it relies 
predominantly on astrometry (plus radial velocities). However, there are two 
major limitations, one practical and one fundamental. The practical problem 
is that applying the method requires high-quality kinematic data in all three 
dimensions (i.e., including parallax and radial velocity). Attempts have been 
made to use only proper motions because they can be very precise, but the 
simulations of iBrown. Dekker fc de Zeeuw ri997l ^ show that doing so leads to 
consistently underestimated ages and overestimates of the initial sizes of the star- 
forming regions. Also, the re sultant kinematic ages are nearly a lways significantly 



less than evolutionary ages (jBrown. Dekker Sz de Zeeuw 19971 ). for reasons that 
are not well understood. 

Even with full 3-D kinematic data, good accuracy is required to get useful re- 
sults and such data are still not available for many groups of interest because they 
are too distant. However, there are a nuni ber of groups of very young stars near 
to the Sun (see Zuckerman Song (20o3 )) that have members with Hi p parcos 



astr ometry. Kinematic ages for these gr oups are summarized in iMakarov f2nn7l ^ 
and [Fernandez. Figueras &: Torra In most cases the kinematic ages are 



greater than ages from pre-main sequence isochrones, but, as discussed below 
(Sec. 14. 2p . the isochrones for PMS stars have significant uncertainties. The 
errors in the measured input quantities are such that going back further than 
~ 20 — 30 Myr is problematic with the kinematic technique. Another limitatio n 
is that some yo ung groups appear to have little net expansion (iMamajek 



Makarov f;2007l );i. so that a time when the members were closest together is ill- 



defined. 

The fundamental limitation on kinematic ages is that stars and clusters en- 
counter massive objects as they orbit in our Galaxy. These objects exert forces 
that disrupt clusters and stellar orbits, creating field stars and leadin g to "disk 



heating," which is the observed increase of vertica l scale height with age (IWielen 19771) . [ 



The time-scale for these encounters is ~ 200 Myr ( Janes Phelps 1994 : Janes. Tillev &: Lvnga 19881 ) .[ 



roughly the Galactic rotation period. Thus one can be reasonably confident that 
the Galactic motions of young stars are undisturbed if they are < 10^ Myr old, 
but for older stars Galactic orbits cannot be reliably traced back in time. 
3.2.1 SUMMARY OF KINEMATIC AGES 
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+ Kinematic age estimation uses few assumptions and involves no stellar mod- 
eling. The only requirements are good three-dimensional space motions, 
meaning high-quality astrometry and radial velocities. 

+ The method is applicable to very young ensembles of stars for which few or 
no other age estimators work well. 

+ The kinematic method can be applied to stars in any mass range. 

— Kinematic traceback can only be done on an ensemble of objects, not individ- 

ual stars. 

— The method works only for young groups, up to r ~ 20 Myr. 

— Accurate parallaxes in particular are required, which are often unavailable for 

faint, low-mass stars of interest. 

3.3 The Lithium Depletion Boundary 

LDB ages have been determined for only a very small number of nearby clusters 
that have been age-dated in other ways. The significance of LDB ages is for the 
overall OC age scale and so a discussion of the lithium depletion boundary is 
given in Section [9. II in the presentation on OC ages. 



4 MODEL-DEPENDENT AGES FOR INDIVIDUAL STARS 

For the most part, the stellar models used to estimate the ages of individual stars 
are the same ones that go into determining OC ages, but there are different ways 
those models are applied because less information is available for the single star. 
OC ages and how they are used to calibrate star ages are discussed in Section [9l 



4.1 Isochrone Placement for Main Sequence Stars 

The principle of deriving an age from placing a star on model isochrones in the 
HRD is straightforward, but in practice it is difficult and problematic. One 
starts with the same isochrones used to fit to cluster CMDs, based on well-tested 
stellar models. However, in fitting a cluster one has not only more stars to 
improve a fit statistically, those stars are also distributed in mass, and it is the 
full behavior of the models over that mass range that is being fit, leading to a 
generally well-constrained situation, depending, of course, on the accuracy of the 
input quantities and the completeness of the physical processes that are included. 

For single stars, that distributed fit is not possible. Moreover, isochrone shapes 
are complex, particularly near the MSTO, and multiple isochrones can pass 
through a given point. Some of that degeneracy is removed if the abundances 
in the star are known, and additional constraints exist if the star is in a binary 
system. The contours and degeneracy of isochrones means that the location of a 
star in a CMD cannot be simply inverted to yield an age. They also mean that 
extracting a most-probable age and realistic uncertainties require special analy- 
sis. This can possibly result in a well-defined age, but sometimes only a limit can 
be established or there is little constraint at all. 

The method starts with best es timates of temperature, i netallicity, and lumi- 
nosity, together with un certai nties. Edvardsson et al. fl993 ]^ Ng fc Bertelli f 19981 ) ,| 



Nordstrom et al. (20041 ) , and Holmberg. Nordstrom Sz Andersen (20091 ) have a^ 



plied this method to the Geneva-Copenhagen survey of ~ 14, 000 nearby FGK 



David R. Soderblom 



14 



stars. One can interpolate among the model isochrones to the location of the 
star, but this introduces biases both because of the complex shapes of isochrones 
and their non-uniform spacing. In particular, input uncertainties in Tefj and Lboi 
may be gaussian, but the uneven spacing of isochrones near the ZAMS (with an 
accelerating increase in Lboi as a star ages) leads to equally probably young- and 
old ages from naive placement of the star, when in fact most stars should fall 
in the dense ZAMS region; this effect biases results toward greater ages and is 
essentially the same as the Malmquist or Eddington bias seen in other areas of 
astronomy. This bias in resultant ages is exacerbated by the fact that stars on 
and near the ZAMS often have indeterminate ages and the more evolved stars 
have well-defined ages (see below). This can also bias apparent ages to higher 
values. Thi s bias can be seen in one of the earliest efforts at precision isochrone 
placement ( Perrin et al. 19771 ). where most derived ages of disk stars exceeded 10 
Gyr and ranged up to ~ 25 Gyr. 

Some recent st udies have continued to estim ate ages from direct isochrone in- 
terpolation (e.g.. i Kartag. Bilir &: Schuster (200 5)). but more sop histicated tech- 
niques have been applied to attempt to overcome these difficulties. IPont fc Ever f2004l ).| 
for example, apply bayesian methods to the F and early-G dwarfs in the Geneva- 
Copenhagen sample, in particular stars with r ~3-15 Gyr, which is to say stars 
that have evolved for a,t least ^ 1/3 of the ir MS lifetime so that they have moved 
away from the ZAMS. IPont k Ever (2004 ) were a ble to lessen biases in the ages , 
although indi vidual ages had lar ge uncertainties. |j0rgensen fc Lindegren (200^ ) 
expanded on iPont &: Ever (2004 ) by starting with a large sample of synthetic 
stars with assigned masses and ages. After adding errors they tested their ability 
to recover the starting ages. They tried several techniques and found bayesian 
methods gave the most accurate net answer with the least error. That uncer- 
tainty was about 0.20 dex, or 50%, on average. Relatively massive and evolved 
stars could be dated to within about 20%. Note that their sample did not go 
as low in mass as 1 Mq, nor did they treat ZAMS stars. Their interest was in 
improving fits to the MSTO in sparse older clusters and they were able to do so 
fo r with remarkably s mall fitting errors. 

Takeda et al. r2007l ) also developed a bayesian method that they applied to a 



large sample of G dwarfs with masses from abo ut .8 to 1.3 Mp), significantly lower 
masses than what either iPont &: Ever (20041') or Ijorgensen &z Lindegren (200,51 ) 
considered. Figured! taken from lTakeda et al. (20071 ). illustrates important as- 
pects of the technique in that well-defined ages result in only some circumstances, 
while in other cases only limits to age can be determined. The net effect of this 
unavoidable effect, as noted, is to add a bias to the resulting age distributions, 
favoring older ages over younger, and more massive stars over the less massive. 

Derived uncertainties in ages from isochrone placement are signifi cant. In their 
analysis of their synthetic dataset, |j0rgensen fc Lindegren (20051 1 were able to 
recover ages to better than 20% only for relatively massive and well-evolved stars. 
Less-evolved stars had ages good to about 50%, and stars near the ZAMS were 
in determinate. 

Takeda et al. (20071 1 got similar uncertain ties, with more ma ssive stars (~ 
1.4Af0) having ages good to ~ 10%. In the iTakeda et al. (20071 1 study, which 
was able to use some of the highest quality data available for FGK stars, only 
2/3 of the stars ended up with a well-defined age in the sense that a clear peak 
was present above the 60% probability level. Of those stars with well-defined 
ages, only 25% had formal errors of less than 1 Gyr, and those are fitting errors 



Fig. [1} Probability dis- 
tribution functions of 
age from Takeda et al. 
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that do not take into account systematic effects in the models. 

The age errors are dominated by uncertainty in T^e, usually taken to be 80- 
100 K. Given a good parallax, the only factor adding to uncertainty in Ly^^i is 
the bolometric correction, and BCs are well-established, at least for stars similar 
to the Sun. But the stellar temperature scale remains fundamentally uncertain, 
in part because we try to characterize what are manifestly non-uniform stellar 
atmospheres with a single parameter that describes the total energy output (Tcs) 
but which does not necessarily describe the temperature structure of the star's 
atmosphere and the conditions under which spectrum lines are formed. Better 
ways of describing the atmospheres of late-type stars are needed in order to better 
relate observations to models. 

4.1.1 SUMMARY OF ISOCHRONE PLACEMENT FOR MS STARS 

-|- Isochrone placement uses models that are based on well-understood stellar 
physics, especially for stars similar to the Sun. 

-I- The method is well-suited to older or more massive stars where other tech- 
niques (such as the empirical methods) are of little use. 

— Not all stars end up with an isochrone placement age even with good quality 

observational data. This effect can bias interpretations of the ages of a 
group of stars. 

— Unresolved binaries will appear to be more luminous (i.e., older) than they 

really are. 

— Error analysis is complex and non-gaussian. 

— Uncertainties are significant even when a well-defined age results, often being 

20-50% before systematic effects are taken into account. 

— Isochrone placement is limited to stars at least as old as ~ 1/3 of their overall 

MS lifetime so that they have moved away from the ZAMS in the CMD. 
This favors more massive stars. 

— Ages and uncertainties are limited significantly by errors in Tgff. Incomplete 

knowledge of composition is also a limitation but generally a lesser one 
if a spectroscopic metallicity is available. The star's helium abundance is 
crucial, of course, but with little information to work with He is assumed 
to be solar. 

— Many different isochroncs can pass through a given location of the HRD. 

Resolving this degeneracy requires at least good knowledge of the star's 
abundances. 

4.2 Isochrone Placement for Pre-Main Sequence Stcirs 

At the most basic level, for PMS stars it is very difficult to be sure one has 
determined an accurate L^oi and T^s, both of which are essential for placing a 
star in an HRD. Star- forming regions remain rich in the residual material from the 
formation process, leaving with it ambiguity about what reddening law to apply, 
how much overall extinction may exist, and how much local extinction may play 
a role, especially since the material closest to stars is clearly not isotropically 
distributed. In some cases it can be difficult to isolate the light of a given star 
from surrounding material or other objects in close proximity. 
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Placing PMS stars on isochrones is generally done for groups of stars that are 
found in a given star-forming region or in an ensemble. Nevertheless, each star 
must be treated individually because of the many factors that affect its location 
in the HRD. Reddening and extinction affect both colors and magnitudes, and 
both can vary significantly in a star-forming region because of the residual mate- 
rial. PMS stars are inherently variable too. Even after correcting for reddening, 
the colors of PMS stars are not reliable temperature indicators. Some PMS stars 
exhibit infrared excesses, indicative of circumstellar material. Others show blue 
excesses and/or excess continuum emission ("veiling") that is attributed to ac- 
cretion. Teff values for PMS stars are best assigned from spectral types and that 
requires additional observations, and the result are imprecise. The relation be- 
tween spectral type and Tefj is likely different for PMS stars than for dwarfs, 
in part due to gravity effects (jLuhman et al. 20031 ) and also because PMS stars 
appear to have highly inhomogeneous atmospheres not well characterized by a 
single temperature. 

Even if one is able to work through these problems to some satisfactory level, 
PMS ages are determined from comparing the locations of stars in an HRD to 
calculated isochrones and mass tracks. These models differ substantially, mainly 
because of differing treatments of convection and also the equations of state and 
the radiative opacities that are used. Uncertainties in the ages of factors of 2-3 
are likely. Figure [2] shows examples of calculated PMS isochrones for low-mass 
stars, with the differences growing in going to the cooler K and M stars which 
form the dominant portion of the T Tauri population. 

The classic c ase study of a sta r-forming region is that of the Orion Nebula 
Cluster fONCV lmilenbrand ^19971 ^ presented a thorough analysis of the available 
photometry and spectroscopy that has recently be e n updated with HST photom- 
etry for many more stars by iDa Rio et al. r2009l ). IPa Rio et al. r2009l ) present 
an analysis that statistically takes into account reddening, accretion effects, and 
so on, plus average luminosities for well-observed stars, yet there remains a sig- 
nificant spread in Lboi at a given Teg. These spreads are typically seen in all 
star- forming regions and are not fully understood. 

A good test of isochrone placement is offered by bina ry systems when consistent 
data of good quality are available for both components. Kraus fc: Hillenbrand (2009bl )| 
found that binaries in Taurus-Auriga, when placed on PMS isochrones, differed 
significantly less in apparent age than did randomly-selected pairs of stars from 
that star-forming region. This was interpreted to mean that stars within binaries 
form with less delay relative to one another compared to the overall time-scale 
for formation within Tau-Aur. Also, the binary components had relative Teg 
and Lhni values co n sistent with the isoch rones used, a blend of models fro m 
Baraffe et al. f 19981 ). Ichabrier et al. f2000l ) and lO'Antona MazzitelU f 19971 ). 



Fig. [2} Examples of 
PMS isochrones. 



Fig. [S HRD for the 
Orion Nebula Cluster. 



However, one can still hope to determine if significant age spreads exist within 
star-forming regions; see Sec. I4.2.2I 

Measured masses from PMS binary star orbits are rare and imprecise, yet are 
essential for calibrating the models. Even the nearest star-forming regions are 
too far away to have trigonometrically-determined distances, although in a few 
cases (Orion and T Tauri itself, notably), radio interferometry has been able to 
greatly improve on that vital quantity. 

As stars approach the ZAMS, their evolution slows and the isochrones crowd 
close together, too close relative to error bars we now have on luminosities, making 
age estimation difficult even if the models are accepted at face value. A few en- 
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sembles of co-moving young stars have been identified (jZuckerman &: Song 20041 ) 



and have members with trigonometric parallaxes from Hipparcos. The ages of 
these ensembles have been estimated from isochrone fitting which works, as for 
star-f orming regions, bec ause a number of stars are present. In a very few cases 
(e.g., iBoden et al. ^200,^ )) it is possible to estimate the age of a single object, in 



part because the "object" is in fact a multiple system, but even then the age is 
imprecise. 

Systematic differences between sets of model isochrones for PMS stars range 
from 20% for the hotter stars (F and G stars, with better-understood physics) 
to as much as factors of 4-5 for late-K and M stars. Aside from that, observa- 
tional uncertainty dominates in comparing stars within star-forming regions (to 
determine age spreads) or between different regions. 

4.2.1 SUMMARY OF ISOCHRONE PLACEMENT FOR PMS STARS 

-|- For PMS stars, isochrone placement uses a property of the stars that is chang- 
ing rapidly and systematically, the luminosity. However, the luminosity 
may not be changing monotonically, given the many factors that influence 
it. 

+ The method is potentially well-suited to more massive stars because differences 
in the models are less. 

— Differences in the available models can be very large, leading to significant 

systematic uncertainty. 

— Placing individual stars in the HRD requires estimates of T^s and I/boi, both 

of which are affected by non-uniform reddening and extinction, and by 
circumstellar material or accretion. 

— Unresolved binaries will appear to be more luminous (i.e., younger) than they 

really are and will add to an apparent age spread. 

— PMS models depend on assumptions of the accretion histories of objects that 

are difficult to verify. 

— Error analysis is difficult because of many poorly-understood factors that lead 

from the observations to the modeled quantities. 

— Isochrone placement for PMS stars is less certain as stars approach the ZAMS 

because isochrones get closer together. 

— Ages and uncertainties are limited significantly by errors in T^s- 

4.2.2 AGE SPREADS IN STAR-FORMING REGIONS Even if the 
absolute age of a star-forming region such as the Orion Nebular Cluster may be 
difficult to establish, we would like to be able to determine the extent, if any, of 
age spreads within such groups. This amounts to having a reliable age ordering 
scheme. The usual method is to examine the spread in L^oi and to compare 
that to expectations from known error sources to see if a residual spread remains 
that is attributable to an age spread. The Orion Nebula Cluster (ONC) and the 
region around it presents a favorable case observationally in that the ONC is 
among the youngest of nearby star-forming regions, helpin g: to make the effect s 
age differences relatively larger. The most recent analysis ( Da Rio et al. 20091 ). 



like those before it, finds no convincing evidence for a spread in any one region 
(see Fig. [3|). There is, however, a progression in age sp atially such that th e 
outer regions of the ONC are older than the central regions ( Da Rio et al. 20091 ). 
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Jeffries et al. (200i( ) studied NGC 2169, a 9 Myr old cluster, and found At < 2.5 



Myr. 

There must be a finite duration over which star formation occurs in a given 
region, the question is measuring it in the presence of many other effects. One 
imagines, though, that binary systems would have very little lag in the formation 
of the two stars and so that binaries in star-forming regions might offer a means 
of assessing the apparent spread within a population with At ~ 0. That ex- 
pectation is borne out for the most part ( Kraus &: Hillenbrand 20()9bl ). although 



inevitably some binaries are conspicuous exceptions for reasons that should help 
u s understand the H RDs of very young groups of stars. 

Palla et al. (20051 ) have applied another means of estimating the age spread 



within the ONC: differences in lithium abundances. Li in PMS stars presents its 
own problems (s ee Sec. 15.31). but a gain differences should be able to be established 



more reliably. (jPalla et al. 2005l l find several low-mass (~ 0.4 Mq) stars that 



appear to be genuine ONC members yet have depleted enough Li to be about 10 
Myr old, an age much greater than the rest of the cluster (~ 1 — 3 Myr). 

Binary systems als o can provide consistency ch ecks on methods of estimating 
ages and age spreads. White &: Hillenbrand (2005! ) show that a probable member 



of the Taurus- Auriga associat ion lacks Li in both compo nents, suggesting an age 
inconsistent with isochrones. Hartigan &: Kenvon (20031 ) also looked at binaries 



in Taurus-Auriga and found that primaries appeared to be systematically older 
than their secondaries for three different sets of PMS evolutionary tracks and 
that some pairs were unusually discrepant, perhaps because of underestimated 
veiling. 

The accretion histories of stars will influence their apparent ages and app arent 

age differences as well, and this has b een addressed b vlFroebrich et al. f2006l) and 

Baraffe. Chabrier &: Gallardo (20091 ) . In particular, Bar affe. Chabrier &: Gallardo (2009l )| 



find that allowing episodic accretion can produce enough spread in Lboi at the 
age of the ONC (~ 1 Myr old) to show an apparent spread of ~ 10 Myr. This 
episodic accret ion also affects stellar radii and, possibly, the existence of the stel- 
lar birthhne of lStahler (19881 ). 

Testing for age spreads within star-forming regions is an instance of determining 
relative ages. iHillenbrand. Bauermeister &: White (20081 1 present a review of this 
subject. 

4.3 Asteroseismology 

Stellar pulsations have long been a means of testing models of stellar physics 
because pulsations probe a star's interior. Isochrone fitting works with the same 
models over the full extent of the HRD but has only the simple surface descriptors 
- Lboi and Tefj - to work with. By contrast, helioseismology has provided such 
a wealth of detailed information and precise physical constraints on the solar 
interior that it has led to fundamental ramifications for particle physics. 

Nearly all types of stars pulsate to some extent; here I concentrate on solar- 
type stars because the oscillations provide a way of constraining a star's age, 
and I focus on just the parts relevant to ages since a steroseismology is it self a 
deep subject worthy of a penetrating examination. See Cunha et al. (20071 ) for a 



recent review of asteroseismology and the physics underlying it. 

For the Sun, ~ 10^ modes have been detected, but we cannot resolve the 
surfaces of stars and so are limited to modes of very low degree. However, those 
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low-order modes probe the deep interior of the star and so set valuable constraints 
on important properties, including, indirectly, the age. The observations needed 
to detect oscillations in stars like the Sun are extraordinarily de manding. The 
Sun's ampli tude is ^ 1 m s~^ in velocity (IKjeldsen et al. 2008al ) and 4 x 10~^ 
in intensity ( Lebreton fc Montalban 20091 ). and these scale roughly in proportion 
to Lboi and inversely with mass. This degree of photometric precision is only 
achieved from space, but ground-based radial velocity observations have yielded 
results for a small number of the brightest FGK dwarfs. 

For the Sun, the primary period seen is ~ 5 min, but the power spectrum shows 
a pattern, with a large frequency spacing Az^ ~ 100//Hz and a small frequency 
spacing 61^ ~ lO^Hz. Av scales as the mean density of the star and is mainly 
sensitive to stellar mass, di^ is sensitive to the gradient of the sound speed which 
reflects the star's evolution. A plot of bv j Av versus IS.v has been shown to be 
a good diagnostic of mass and age for solar-type stars that should lead to ag e 
accuracies of ^ 1 % (Oti Flo ranes. Christensen-Dalsgaard fc Thomp son (20051 ): 
Mazumdar (20051 ): iKieldsen. B edding h Christensen-Dalsgaard (20081^)). Detect-| 



ing these frequency spacings requires a near-continuous data stream for about a 
week, meaning, when done from the ground, a multi-site campaign to eliminate 
diurnal aliases in the power spectrum. At the same time, the signal sought is by 
its nature periodic and the mode lifetimes are long so that summing observations 
can yield the needed signal-to-noise. 

Ages from asteroseismology use the same models applied to isochrone place- 
ment and so have the same model vulnerabilities. As for isochrone placement, 
it is important to constrain basic stellar parameters as well as possible, par- 
ticularly the composition. In some cases it may be possible to independently 
measure a star's radius (via interferometry) , mass (from a binary orbit), or 
mean density (from a planetary transit). The great value of asteroseismology 
is that it yields more significant constraints on the models, and can include 
the He abundance and the s ize of the convective envelope or convective core 
(|Lebreton fc Montalban 20091 ). 

Because the observations are so resource-intensive only a small number of stars 
have been obse rved (I am aware of five with asteroseismological ages; see, e.g., 
Vauclair (20091 )). In addition, a full fit to the observations for any one star 



requires a custom-built stellar model tuned to each particular case. As discussed 
below (Sec. [TT]l . this is about to change significantly because of the CoRoT and, 
especially, Kepler missions, which will produce high-quality asteroseismological 
signatures for thousands of stars. What has been a boutique business is about to 
go past retail to wholesale. Among the few stars with detected seismic frequencies, 
even fewer have had ages determined. However, the ages can be very precise, and, 
since they use information from much of the star, are likely to be more accurate 
as well. 

With a sufficiently long time-span (months to a year), asteroseismology is ca- 
pable of detecting aspects of the star's interior that arise from convection, which 
is to say the depth of the convective envelope. Only Kepler is likely to yield such 
data, but it can provide for the first time critical parameters of interior structure 
for convective stars other than the Sun. 

4.3.1 SUMMARY OF AGES FROM ASTEROSEISMOLOGY 



-|- Asteroseismology provides physical constraints for the star observed that are 
available in no other way. 
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+ Asteroseismology is especially effective for older MS stars because the increase 
of He in the core increases the density, leading to an increase in the sound 
speed, which in reflected in the oscillation frequencies. It is these stars for 
which the ages can be so intractable otherwise. 

+ The ages that result from a full modeling effort are probably accurate to 
~ 10%, significantly better than any other age-dating method. 

+ Asteroseismology can be applied to individual stars at any evolutionary state. 

— Asteroseismology uses the same stellar models used for isochrone placement 

(for individual stars) and isochrone fitting (for clusters) and has the same 
limitations. At the same time, the stars observed are very similar to the Sun 
and so the models and the underlying physics should be sound. Efforts are 
underway to better u nderstand the limitations and problems in the models 
(|Lebreton et al. 200^ ). 

— Obtaining the needed observations for a particular field star is extremely re- 

source intensive, requiring the largest telescopes and best spectrographs in 
a coordinated campaign, followed by modeling that is tailored to the source. 

— Our understanding of the physics of PMS and low-mass stars is not yet suffi- 

cient to predict or analyze their oscillation properties. 

5 EMPIRICAL AGE RELATIONS 

Empirical age relations provide an observed change of a quantity with age that can 
be calibrated with clusters or stars that have ages determined in other ways. Also, 
a plausible physical scenario exists to provide a framework for the observations, 
even if we do not fully understand the physical processes and cannot predict 
the behavior of the property. Because empirical indicators are calibrated against 
objects with model-dependent ages (clusters mostly), they are secondary to the 
model-dependent methods. 

The empirical relations discussed here apply to low-mass stars (main sequence 
and PMS) and this is because they all have their origin in the presence of a con- 
vective zone (CZ) in such stars. In such stars, convection and rotation (especially 
latitudinal differential rotation) interact to create complex motions, and the ma- 
terial in the CZ is ionized and therefore conducting. These motions can then 
regenerate a seed magnetic field to reinforce and amplify it: this is the dynamo 
mechanism. Going further, we observe the Sun to have an ionized wind, and 
that enables the magnetic field to force corotation of the wind to well beyond the 
Sun's surface, leading to angular momentum (AM) loss. Thus in this picture the 
gradual spindown of a star like the Sun is inevitable. Also, stars that start with 
higher-than-average rotation rates produce stronger magnetic fields and so lose 
AM faster. In other words, there is a feedback mechanism that leads to an initial 
spread in rotation rates among a coeval population converging over time. 

The magnetic fields that are generated through the dynamo mechanism mani- 
fest themselves as activity in a variety of forms. Non-thermal energy is deposited 
into the chromosphere and corona. Properties of a star's chromosphere can be 
observed by obtaining spectra in the cores of very strong absorption features such 
as Ha, the Mg ii h and k lines, the Ca ii H and K lines, or the Ca ii infrared 
triplet. At those wavelengths the stellar photosphere is suppressed by the ab- 
sorption and the weak chromosphere can be seen. Stellar coronae are seen as soft 
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x-ray emission. 

This scenario forms the paradigm that guides the study of spindown, the decHne 
of activity, and the relation between the two in late-type stars. It is important 
to emphasize that this picture is an hypothesis that has been partially tested but 
not completely. The processes that feed energy into the chromosphere and corona 
are not understood even for the Sun. However, this scenario is consistent with the 
observations available so far (but see a comment below). Another consequence 
of the CZ is, we believe, the depletion of a star's surface Li abundance over time 
and that provides an additional empirical age indicator. The spots and magnetic 
activity that are associated with these processes lead to photometric variability 
as well, and that is briefly described. 

This subj ect o f empirical age ind icators has its own history, but th e papers by 



inis sub.i ect o r empirical age ind icators nas its own nistory, but tn e papers by 

Kraft r 19671 ) and lSkumanich (19721 ) have inspired much that followed. ISkumanich (197^ ) 



used rotation, activity, and Li data then available for the few nearest OCs and 
for the Sun and noted that rotation and activity both appeared to decline in pro- 
portion to T~^/^, while Li depletion followed an exponential. The rotation data 
available at the time was from vs'mi measurements, but in recent years precise 
photometry has yielded an abundance of accurate rotational periods for solar- 
type stars, especially in younger clusters, and this has enabled a much better 
view of the phenomena taking place. 

5.1 Rotational Spindown and Gyrochronology 

The decline of rotation with age for solar-type stars is well-established from obser- 
vations of clusters and field stars. Because of the feedback scenario just outlined, 
this spindown is thought to be deterministic, not probabilistic, and there is a con- 
vergence in rotation rates with time among coeval populations. Still, the details 
of AM loss remain to be investigated and there are poorly-understood aspects 
of spindown that challenge our picture, particularly for PMS and ZAMS stars. 
Spindo wn and AM evolution have been addre ssed i n many reviews; recent ones 
in clude lBarnes/2nn9l '). llrwin fc Bouvier (20Q(h . and lMamajek f 20091 ). 



Barnes (2Qm has calibrated a relation betwe en rotation period, color, and age , 



using data for clusters, and this was updated by Mamajek Sz Hillenbrand (20081 ). 



These calibrations involve some assumptions that affect how they can be applied 
to individual stars. 

1. As noted, OCs older than ~ 500 Myr are relatively rare and so tend to be 
distant. This inhibits testing the models in the realm they are most needed: 
among the older stars where one wishes to estimate ages of individual field 
stars. This can be addressed to some extent by finding nearby binaries 
with intermediate-mass primaries and solar-type secondaries so that the 
system's age can be established from isochrone placement, but there are 
not very many such systems to work with. 

2. Barnes (20071 ) characterizes each cluster with a mean relation fitted through 



the stars with the longest periods. Indeed, most members in young clus- 
ters are found on his "I" sequence, but in any young cluster there is a 
substantial fraction of more rapidly-rotating stars (his "C" sequence), in 
some cases spinning at 100 times the rate of the slower stars (Fig. 
These ultra-fast rotators (UFRs) are themselves a fascinating phenomenon 
with important implications for AM partition and evolution among PMS 
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stars. They are bona fide cluster members and add substantial cosmic 
"noise" and uncertaint y to the calibration. This can be seen clearly in 



Irwin &: Bouvier who illustrate the available cluster data. 



3. In addition to adding uncertainty to age determinations, the scatter seen 
in young clusters means that ages estimated from Pmt for young stars 
will have a bias. The gyroch ronology calibrations of Barnes (200?! ) and 
Mamaiek fc Hillenbrand (20081 ) are based on means for the slower rotators 
in the various clusters. Those same relations, if applied to one of those 
clusters, would produce a distribution of ages that is asymmetric with a 
pronounced tail toward stars that appear to be younger than most of the 
cluster; these are the UFRs. For an individual star one does not know 
where in such a distribution the star's AM lies. 

4. This spread in AM appears to have converged by an age of ~ 500 Myr, as 
shown by observations of rota tion periods in F and G sta rs in the Hyades 
and Coma Berenices clusters ( Collier Cameron et al. 20091 ). 

5. Rotation rates at a given age are highly mass-sensitive, particularly above 

I Mq on the main sequence (PMS rates tend to be more independent of 
mass; see NGC 2264 in Fig. U]). Calibrations are generally done in terms 
of stellar color as a proxy for mass, but the color-mass relation is different 
for different compositions and this adds systematic uncertainty. 

6. Obtaining the P^ot distribution in clusters is relatively straightforward in 
that young stars tend to show significant photometric modulation from 
surface inhomogeneities (spots) and clusters cover modest regions of the 
sky, making it possible to get the needed observations over a month or so. 
By comparison, the older stars of the field must be observed individually 
and the amplitude of variation tends to be much less. Rotation periods in 
young stars range from a few hours to ~ 2 weeks, while older stars like the 
Sun have periods of a month or more. 

7. Our knowledge of rotation from Pj-ot is also biased because only some stars 
show periodic variability at a given time. They are generally the more 
rapidly rotating and more active stars, but even among those not all stars 
show rotational modulation all the time. A higher-contrast signal (com- 
pared to broad-band measures) can be seen by using observations of a 
chromospheric feature such as Ca ii H and K, although at the cost of signal- 
to-noise. Even then luck plays a role, as illustrated by observations of Mg 

II h and k in the Sun with the SORCE mission: For the most part the 
rotational modulation of the hk emission was unmistakable, but at times 



9. 



little or no modulation could be seen (jSnow et al. 20051 ). 

Rotation periods, when available, are clearly to be preferred as a precise 
and aspect-independent measure of rotation, but use of v sin i data is rea- 
sonable as long as rotation dominates the line broadening. In such cases 
the errors are modest, and the statistics of rotation is such that most stars 
show a vs'mi that is at or near the true e quatorial velocity. For ex ample, 
the distribution of wsini in the Pleiades ( Soderblom et al. 1993bl ). when 
converted to Prot, is just like that from Pr„t fo r the Pleiades or for other 
clusters of the same age (llrwin fc Bouvier 20091 ). 

As with any astronomical phenomenon, we can only observe the OCs that 
nature makes available to us. The available data is impressive and growing 



Fig. m Spread in P^ot 
with mass for represen- 
tative clusters. 
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(jlrwin k Bouvier 20091 ^. yet it is still difficult to assess how typical or not 



an individual cluster may be. An even broader sampling of the nearby 
clusters is needed. 

10. Using rotation to derive an age rests on the critical assumption that AM 
loss and rotational convergence are inevitable and that assumption remains 
untested. It is called into question by the recent finding of a star in M67 (age 
4 Gyr) that rotates at about twice the rate of other stars of the same color , 
which rotate at the solar rate, as expected ( Reiners &: Giampapa 20091 ). 



This observation needs to be confirmed (the object could be a binary, for 
instance), but an unambiguous spread in rotation at such a late age would 
ruin the carefully constructed picture we have created. 



The rotation-age calibrations of lBarnes (2007l 'l and lMamaiek &: Hillenbrand f2008l H 



distill knowledge from clusters into mean relations. It is, of course, helpful to take 
advantage of data for many stars in a cluster to define a mean relation, but the 
scatter seen indicates how much different in age a star could seem to be with- 
out our realizing it. The scatt er seen is not statistical and not an error; it is 
real and potentially confusing. Collier Cameron et al. (20091 ) show an example 



of how much spread in apparent age can exist in a coeval population: For the 
well-behaved Hyades and Coma clusters it is ±40%. In other words, there is an 
inherent liability in applying knowledge of a group phenomenon to individuals. 
At the same time, it should be possible to use these calibrations to obtain the 
age of a young ensemble of, say, 10 or more stars. 
5.1.1 SUMMARY OF GYROCHRONOLOGY 

+ Rotation periods can generally be measured with good precision and accuracy. 
Also, rotation is not subject to variability except for the relatively minor 
effect of differential rotation. 

+ The dominant fraction of post-ZAMS and older stars (r > 500 Myr) are 
reasonably well behaved and exhibit a well-defined trend of rotation with 
age. 

+ A plausible and partially-tested scenario exists to explain the observed spin- 
down and convergence in rotation with age. 

— No age can be estimated if the rotation period is not detected, and that limits 

gyrochronology to solar-type stars (because they have spots) and favors 
younger stars for the same reason. Late-F and early-G dwarfs tend to 
exhibit less of a rotational signature as well; note the lack of stars at the 
bluer colors in Figured 

— Obtaining rotation periods is observationally intensive and is not always pos- 

sible if the level of photometric variations is too low. Low amplitude can 
occur at all levels of activity and can bias the results. 

— In some cases an alias of the true rotation period may be measured instead of 

the true period. 

— Differential rotation on the stellar surface can mean that the measured rotation 

period is not the equatorial value. However, experience so far shows that 
the uncertainty this introduces is modest compared to other effects. 

— Significant scatter in Prot exists among stars in a cl uster which ar e of the same 

age, particularly for r < 500 Myr (Fig. 2]). iBarnes (20071 ) argues that 
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the stars in clusters that are on his "I" sequence follow the r^^/^ law that 
appears to hold for older objects, but it is not possible to tell if a given 
individual star meets that criterion or not. Given the way in which the 
rotation-age relation is calibrated (see above), the resulting ages can have a 
bias because the rapid rotators that are present at young ages are predicted 
to be much younger than they really are. 

Different clusters of the s ame or similar ages do not show identical distribu- 
tions of Prot with mass (jlrwin fc Bouvier 20091 '). adding additional inherent 



uncertainty. There may be other parameters (such as composition) that 
current calibrations do not take into account, or initial conditions may 
matter even after stars reach the ZAMS. 

5.1.2 RECOMMENDATION Despite these known and potential flaws, 
using Pj-ot to estimate an age for a G or K dwarf is one of the better-calibrated and 
better-understood methods that is available. Paradigms are our creations to guide 
our understanding - they do not necessarily represent reality - but the rotation- 
age paradigm has withstood well the tests of it that have been made, given the 
admonitions discussed above. The calibration of Mamaiek &: Hillenbrand (20081 ) 



for F, G, and K stars is to be preferred because it appears to avoid a systematic 
error that Barnes included. 



5.2 Activity and Age 

As mentioned, ISkumanich f;i972l ) first suggested that a r power law may 
relate both rotation and activity to age. In particular, he used measures of Ca ii 
H and K chromospheric activity because those data were available, but there are 
other indicators of activity that can be used as well. The Skumanich relations 
also imply that activity is proportional to rotation and that is indeed observed 
but with two critical limitations. First, the activity we see may have its origins 
in rotation, but the chain of processes that lead from spin to magnetic field 
generation to observed activity is a long one with poorly understood links. One 
result of that chain is that activity can be highly variable for a given star. Second, 
for reasons not well understood, activity saturates at high rotation rates. This 
makes qualitative sense in that activity cannot remain linearly proportional to 
rotation indefinitely, but it adds additional uncertainty. 

What is broadly called activity is observed in a number of forms. All of it is non- 
thermal emission from the star and is much weaker in flux than the photosphere. 
As a result, activity is best observed at wavelengths where the photosphere is weak 
or suppressed, such as in x-rays, the ultraviolet, the cores of strong absorption 
lines, and in the radio regime. For some of these wavelengths (x-ray and radio) 
the signature of activity is unambiguous in that the thermal background is non- 
existent. However, the most easily observed indicators of activity are the cores 
of strong lines in optical spectra and in those cases disentangling the flux due 
to activity from the background is not so clear-cut. In addition, it is not always 
clear how to parameterize the activity, with some observers preferring the flux 
in, say, soft x-rays. Here I work with normalized indices, in which the observed 
flux is divided the stellar bolometric flux to derive a flux ratio such as Rx or 
RuK because these indices appear to produce consistent results across a range of 
spectral types. 
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5.2.1 Ca II CHROMOSPHERIC EMISSION By far the most frequently 
observed manifestation of activity is the chromospheric reversals seen in the cores 
of resonance lines in optical spectra, particularly Ca II H and K. Here I will focus 
primarily on FGK stars but HK data exist also for M dwarfs. The use of Ca H and 
K as an activity indicator for stars was pioneered by O lin Wilson, first with photq - 
graphic spectra and then with a spectrophotometer. Vaughan &: Preston (198(3 ) 
carried Wilson's work forward with a dedicated instrument that has been sem- 
inal for understanding the Sun in the context of the stars and vice versa. The 
HK project at Mt. Wilson has led to long-term monitoring of stellar activity 
cycles for nearly a half-century and has included a volume-limited survey as well 
as observations of a few clusters accessible with their telescope (the Hyades pri- 
marily). Many papers and reviews have described the HK observations and their 
reduction to a physically relevant quantity, -R^j^, the ratio of the HK flux to 
Lhni, corrected for p hotospheric light in t he bandpasses (hence the prime; see 
Noves et al. (1984 )). Baliunas et al. (19951 ) provide a summary of the Mt. Wil- 



son program. 

Observing HK in FGK stars has a very significant advantage over attempting to 
detect Prot: A single low-resolution {R ~ 2000) spectrum suffices for estimating 
the mean activity level. In addition, HK emission can be detected even in older 
stars that lack enough contrast to make their rotation detectable from modulation 
of surface inhomogeneities. That makes it possible to measure HK in thousands 
of field stars. However, HK activity is known to vary on a number of time-scales, 
from hours (fiares), to days and weeks (rotation), to years (long-term activity 
cycles). Each of these regimes is itself a subject of study for understanding stellar 
physics and the dynamo mechanism, but here our concern is that variability 
adds uncertainty in defining what the mean activity level of a given star is. 
Still, these variations are of fairly low amplitude compared to the overall drop of 
activity with age. Data for stars observed many times at Mt. Wilson show that 
one or two HK measureme nts produce an H K equivalent width that is within 
8% of the long-term mean (ISoderblom 198,51 ^. This leads to about a 0.08 dex 
chang e in log-R^^r for an older solar-type star and less for cooler or younger 
stars. IWright et al. (2004 ) compared Mt. Wilson observations to a similar index 
derived from high-resolution spectra for the same stars and saw a 13% scatter. 
Multiple, exte nded measurement s of chromospherically-quiet stars show a scatter 
of about 6% ( Wright et al. 2o"o3 ): this represents the fioor to inherent variation 
in the best cases. As another example, Figure P shows long-term observations 
of HK in the Sun ( Hall. Lockwood Sz Skiff 20071 ) obtained at Lowell as part of 
a stellar monitoring program. Short-term variations can be substantial, and the 

solar cycle is easily seen. Fig. [5} HK measures 

Following the initial look at the age d ependence of activitv bv Skumanich (19721 ) .[for the Sun at different 



Soderblom. Duncan Sz Johnson (19911 ) undertook a more thorough calibration, times over a period of 



confirming that the r~^/^ relation appeared to work well for clusters and binary years, 
companions o f more massive stars, with the saturat ion turn-over at young ages 
noted above. Soderblom. Duncan &: Johnson (199ll ) also applied an additional 
method that took advantage of the unbiased (essentially volume-limited) nature 
of the sample that was available. That method used the Sun to anchor the 
relation at 4.5 Gyr and assumed HK emission declines monotonically with age 
(variability is not critical as long as the age trend dominates) and that the Galac- 
tic star formation rate in the vicinity of the Sun has been cons tant. The result, 
after correcting for "disk heating" in the Galaxy ( Wielen 19771 ). was essentially 
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the same as the directly-cahbrated relation, suggesting both methods are reason- 
ably sound and that the chromospheric data are consistent with a constant star 
formation rate. 

Just as for rotation, HK measures in young (< 500 Myr) clusters show lots of 
scatter that is inherent (Fig. [6]). This scatter is related to the spread in rotation Fig. [6l HK in clusters, 
and it exceeds the differences in cluster averages so that there is significant overlap 
in apparent age for stars with the same HK level. As with rotation, this scatter 
is largely absent in the Hyades and it is assumed to reflect the convergence in 
rotation that has taken place. For older stars, the HK emission gets weak and it 
is not clear which portion of the line to attribute to just the chromosphere. 

Those early analyses assumed all the stars were of solar composition and that 
there was no dependence of the HK emission on composition. That assump- 
tion is certainly not true for the oldest stars in the solar neighborhood. Lower 
metallicity may affect the dynamo itself by changing a star's structure, but, at a 
minimum, decreasing [Fe/H] in a star means decreasing [Ca/H] and so changing 
the line formation conditions. The net result is to make a metal-poor star look 
more active and younger tha n it r eally is because the line c ores are not as deep. 
Rocha-Pinto &: Maciel r2000l ) and iRocha-Pinto et al. have calibrated this 



effect. They also showed that high levels of activity may affect Stromgren pho- 
tometry of stars and thus ages determined in that fashion as well. 

In addition to using HK to estimate ages, an essential question has been the 
overall form of the relation and its implications both for ages and for stellar 
physics. A t~^/^ power law may be mathematically simple yet physically incor- 
rect, and the question arises in part because the distribution of HK strengths for 
solar-type stars shows a "Va ughan-Preston gap" that corresponds to an age of 
~l-2 Gyr. Pace et al. f;2009l l have addressed this issue by obtaining HK obser - 
vati ons of some memb ers of clusters in that age range. IPace fc Pasquini f;2004l ) 
and Pace et al. (20091 ) suggest that HK emission ceases to drop after '^1.5 Gyr 
because of their observations of HK in several clusters near 1-2 Gyr age. Verifying 
and calibrating the continuing decline of activity in older stars is difficult, both 
because of the scarcity of old clusters and because the chromospheric emission is 
weak and so high signal-to-noise spectra are required. However, both the Sun and 
the a Centauri system, at 4-5 Gyr, show weaker HK emission than 1 Gyr old clus- 
ters, and M67 - also about 4 Gyr old ~ also shows weak er HK emission (Fig. El) 
but also a spread in HK that exceeds observational error ( Giampapa et al. 20061 ) . 
There is still much we do not know about the evolution of rotation and activity 
in stars like the Sun. 

5.2.2 OTHER CHROMOSPHERIC INDICATORS Other useful mea- 
su res of chromospheric activity in late-type stars include Ha and Mg ii h and 
k. Lvra fc Porto de Mello (20051 ) present observations and an analysis of the age 
dependence of Ha in solar- type stars, an effort started by Herbig (19851 ). Ha is 
not inherently as deep a feature as the H and K lines and so the chromospheric 
component is not as easily detected, requiring both good resolution and high 
signal-to-noise. The Mg ii h and k lines are very similar to the Ca ii H and K 
lines and formed under similar conditions. ICa,rdiui CassatellTTiH^ ) show a 
study of Mg II hk emission versus age. 

Both Ha and Mg ii hk lack the large datasets that are available for Ca II 
HK. However, in some cases those indices may be available for a star of interest 
when HK is lacking. Given that many fewer cluster observations are available to 
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calibrate Ha and Mg ii as a function of age, the best use of them may be to first 
transform them into an effective HK index using mean relations and then use the 
HK versus age calibration. 

5.2.3 CORONAL EMISSION The coronae of stars like the Sun reveal 
themselves in soft x-rays. This can provide a high-contrast signal because stellar 
emission at those wavelengths is entirely non-thermal, yet the signal is a weak 
one and difficult to detect except for nearby and young stars. The significant 
strength of x-ray emission in PMS and ZAMS late-type stars has made survey 
missions such as Rosat a n effective means of ident ifying previously unknown very 



young stars in the field (jZuckerman Sz Song 20041 1 



The limited data available means th at there is not an independent ca libration 
of x-ray emission versus age, although iGiidel. Guinan &: Skinner studied 



a small sample of solar-type stars with ages determined in other ways and found 
Lx oc T~^'^, a much steeper age dependence than the t~^/^ seen for HK emission. 
The situation for x-rays also involves substantial changes in coronal structure 
with time as the temperatures change, and so the measured decline depends 
on the x-ray wavelength s included in the observation. It has also been shown 
(jSterzik fc Schmitt 19971 ) ";hat x-ray emission and Ca HK emission are realted to 



one another with a power l aw over nearly four decades of x- rav flux, with a scatter 
of about 0.06 dex in -Rjjk- iMamajek &: Hillenbrand" !(200i) show a mean relation 
to convert x-ray flux into HK. There is no advantage to using x-ray emission as 
an age indicator if HK data are available, and, in addition, the x-ray fluxes are 
significantly more variable that is HK, adding u ncertainty. Like chrq mospheric 
emission, x-ray emission saturates at high levels ( Pizzolato et al. 20031 ). 



5.2.4 SUMMARY OF USING ACTIVITY TO ESTIMATE AGE 

+ Obtaining an activity index for a star is observationally straightforward and 
so large datasets exist. More activity data exist than for any other age 
indicator, particularly for field stars. 

+ Ca II H and K are easily observed from the ground and only moderate resolu- 
tion is needed to get a measure of the activity. 

+ Although activity is variable on nearly all time-scales, for main sequence stars 
that variability is generally much less than the broad decline of activity with 

age and so even a single observation provides a reasonable estimate of the 

mean level. Ages can be estimated to about 0.2 dex ( Sod erblom. Duncan &: Johnson 199ll ).| 

+ X-ray emission can also be used. It is generally detectable only for young 
stars, but it presents a high-contrast signature. 

+ For a given age, activity is only moderately mass-dependent ( Mamajek &: Hillenbrand 20081 ) .[ 



Activity arises from the interaction of convection and rotation (as far as we 
know) and so it follows from the rotation, but with significant variability 
added. In other words, activity is secondary to rotation, and a given rota- 
tion rate in a star can result in a range of activity. This paradigm needs 
further study and verification, and could be wrong in important respects. 
Also, activity in these forms is only seen in convective stars, about F6V 
and later. 

Young (PMS and ZAMS) stars exhibit a broad range of activity at a given age. 
This prevents activity from being used to date single stars, and ensembles 
of at least 10-20 stars are needed to define a mean relation. 
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Activity can be difficult to detect in older stars, and the portion of the signal 
to attribute to age (as opposed to a basal level expected to be present in 
a minimally active star) is uncertain. We know the Sun has experienced 
periods such as the Maunder minimum in which spots are rare, but we do 
not know the activity level one would observe in those cases. 

Currently-measured activity indices (especially Ca ii H and K) probably con- 
tain a sensitivity to metallicity that has not been fully taken into account. 
One result is that metal-poor stars can appear younger than they really 
are. 

Activity indices such as -Rhk generally measured in a way that does not 
take into ac count rotational broadening of lines , leading to systematic mis- 
estimates. ISchroder. Reiners fc Schmitt show how to account for 



this effect. 

— Activity is inherently variable on all time-scales, and so extended measure- 

ments are needed to get the best estimate of a mean level. 

— Activity saturates at high rotation rates so that activity becomes only weakly 

dependent on Prot and so is a poor age indicator in that regime. 

— Comprehensive data are available only for Ca ii H and K. Other activity 

measures have not been calibrated against age. 

5.2.5 RECOMMENDATION Activity observations need to be first re- 
duced to an index such as -Rhk ^^^^ normalized to the star's luminosity (see 
Wright et al. (20041 ) for the procedure for Ca HK). There are extensive obser- 



vations for Ca II H and K, including stars in clusters, but many fewer data for 
other activity indices as well as fewer observations of clusters. That makes it 
advantageous to first convert an activity observation into the HK equivalent with 

a mean relation such as that shown in Figure 

Many HK-age relations have been published. In Soderblom. Duncan &: Johnson (199ll )| 



an activity-age relation was presented which was calculated for a constant star 
formation rate that incorporated a correction for disk heating. Because this was 
not analytical, Soderblom. Duncan &: Johnson (199ll ) did not provide a formula 



to go with that curve (shown in their Fig. 10), but the following relation is a 
close approximation for logi^^jj- values from —5.1 to —4.3: 

logTHK = -8568.2 - 7037.6X - 2164.6x2 - 295.76x^ - 15.144x^ (1) 

where x = log PjjK ^^'^ ''"hk is in yr. I emphasize that this formula is a convenient 
mathematical fit to the means for clusters and to the Sun; it is not a physical 
relationship a nd it does not take into account the considerable spread seen at 
any one age. I Wright et al. r2004l ) present a simil ar polynomial; the d ifferences 



between the two are < 5% above 2 Gyr, but the Wright et al. f2004l ) ages are 



systematically greater t han Soderblom. Duncan Sz Johnson fl99ll ) below that. 



The turnover shown in ISoderblom. Duncan Sz Johnson fl99ll ) corrects for this 



effect, which is due to saturation of the HK emission. 



5.3 Lithium Depletion 

We know that stars like the Sun deplete their surface Li abundances over time 
because primordial solar system material has about 200 times the Li content of 
the present solar surface. Just as for rotation and activity, young OCs show that 
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young solar-type stars have much more Li than MS stars and also that there is 
substantial and real scatter in Li a t those young ages (Fig . [7]). Some of that 

scatter may correlate with rotation ( Soderblom et al. 1993al ). Fig. [3 Li in clusters to 

Li is observed in stars as absorption at the 6708 A resonance doublet. This show scatter, 
feature is analogous to the Na D lines or Ca ii H and K but is unresolved. The 
weakness of this feature in the Sun {Wx < 2 mA) demonstrates the element's 
paucity in the solar atmosphere. The 6708 A feature is due to Li i. The ion- 
ization potential of Li is low and so most of it is in the form of Li ii which is 
invisible. Because of this, the formation of the 6708 A feature is highly sensi- 
tive to Tpff. Non-LTE effects can also be significant, particularly for Li-rich stars 
(iLind. Asplund k Barklem 2009l l. 

Detecting the Li 6708 A feature in young stars is easy {Wx > 100 mA) and 
indeed the presence of the feature is used as a discriminant for T Tauri stars. A 
secondary Li feature is sometimes seen in Li-rich stars at 6103 A. Detecting the 
Li feature in Li-poor stars like the Sun is very difficult because the 6708 A feature 
is blended with many other lines of comparable strength. Extracting the solar L i 
abundance was a to ur de force, both in the observations (iBrault &: Miiller 19751 ) 
a nd in the analysis (Miiller. Pevtremann &: de la Reza 19751 ). 



Sestito k Randich f2005l ) summarize and analyze the available Li data for OC 



stars. These data are consistent with all the clusters have formed with a uniform 
initial Li abundance of about 3.2 (on a logarithmic scale where H = 12). They 
show that most Li depletion occurs before stars reach the ZAMS, with some 
continuing depletion up to at most 1-2 Gyr. The OC data also show that Li 
depletion is not a continuous process and cannot be described by a simple r~" 
type of law. 

None of the models proposed for Li depletion can fully explain the observa- 
tions, but there are several factors that may contribute to fundamental uncer- 
tainty in interpreting Li abundances. First, because of the temperature sensitivity 
just noted, the observed Wx could be altered by stellar surface inhomogeneities 
(spots). In other words, starspots could lea d to an apparent spread in Li where 
none really exists. ISoderblom et al. (1993al ) show that the degree of spottedness 
would have to be extreme to produce the scatter in Li seen in the Pleiades, for 
instance. Also, T Tauri stars, which we observe to hav e very high levels of spot - 
tedness, show less Li sc atter than ZAMS clusters do ( Sestito Sz Randich 20051 ). 
Soderblom et al. (1993al ;i looked for Li variability that might arise from spots and 



saw none, but such an effect has not been fully ruled out. 

A second potential effect that could add scatter would be late deposition of 
gas-depleted material onto a star's surface after planets form. The apparent Li 
abundance would be especially enhanced because the star's surface Li content 
is low to begin with. This process would be expected to lead to higher Li en- 
hancements in, say, a late-F dwarf compared to a K dwarf because the F star has 
much less mass in its convective envelope. Searches for this effect have seen none 
and there is no conclusive evidence that stellar Li abundances are related to the 
p resence of planets. 

Mentuch et al. (20081 ) used observations of Li in young, nearby associations and 



fitted them to calculated Li isochrones to derive ages. This has the advantage of 
using data from a number of stars to arrive a t a net age but makes the results 
model-dependent, do Nascimento et al. (20091 ) go much further and use models 
of Li depletion with non-standard mixing processes to derive ages for stars as old 
as the Sun. Doing so requires assumptions about the histories of the stars that 
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cannot be tested, particularly their angular momentum evolution. 

The depletion of Li in low-mass stars is a significant subject in itself and cannot 
be fully covered here. The mechanism is believed to be convection (or convec- 
tive overshoot) combined with non-standard mixing and observations of Li are 
important for probing fundamental aspects of solar-type stars. However, the use 
of Li for estimating ages is problematic for a number of reasons in addition to 
those just noted. Late-F and early-G dwarfs deplete little Li and the scatter 
seen within clusters exceeds any age trend. Li may work better for late-G and 
K dwarfs, provided that a group of 10-20 stars is examined to average out the 
inherent scatter. 

5.3.1 SUMMARY OF LITHIUM DEPLETION 

-|- The primary Li feature at 6708 A is easy to detect and measure in young solar- 
type stars, having an equivalent width of 100 mA or more and lying near 
the peak of sensitivity for many modern spectrographs. Good resolution 
(at least R « 20, 000) and good signal-to-noise are needed, however. 

-|- The trend of declining Li abundance with age is strong for PMS stars, partic- 
ularly K dwarfs. 

-|- Many clusters have been observed to provide a calibration. 

-|- Clusters appear to form with a uniform initial Li abundance. 

-|- The youngest stars (T Tauris) all appea r to be rich in Li, confirrn ing our 
expectation. There may be exceptions ( White Sz Hillenbrand 20051 ). but in 



such small numbers that special circumstances may apply. 

-|- Young FGK stars appear to always show a strong Li feature, making its pres- 
ence a discriminator of youth at least, even if a quantitative age is difficult 
to establish. 

— Little change is seen in Li with age for late-F and early-G stars. 

— The physical processes that drive Li depleti on remain poorly underst ood. 

Models can reproduce the solar abundance ( Charbonnel &: Talon 20051 ) but 



detecting Li in stars as old or older than the Sun is very difficult and so 
there are few constraints on the theory at present. 

PMS and ZAMS stars of a common age exhibit a large range of apparent Li 
abundance for reasons that are not understood. As is the case for rotation 
and activity, a single star can be placed among clusters with a broad range 
of age, and so an ensemble of 10-20 stars is needed to define a mean relation. 

Converting an observed Li equivalent width to an abundance is highly temper- 
ature sensitive and can require a significant non-LTE correction. Despite 
this, one can work with just the equivalent width and color to compare an 
ensemble of stars to clusters. 

The apparent Li abundance of a star may reflect unappreciated processes that 
could add Li-rich material to the surface (left over from planet formation) 
or co uld remove the star's surface la yers during an early phase of high mass 



loss (jSackmann Sz Boothrovd 20031 1. Starspots may alter the equivalent 
width of the Li feat ure for a given abun dance. PMS magnetic fields may 
inhibit Li depletion ( Ventura et al. 19981 ). 
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5.4 Photometric Variability 

All stars vary to some extent. Stars like the Sun are known to vary on many 
time-scales, and the dominant cause is surface inhomogeneities (spots). Younger 
stars rotate faster than older stars and have higher activity levels and greater 
spot coverage. As a result, there is an overall trend of the degree of photometric 
variability with age. This is not yet well enough characterized to be useful as an 
age indicator, and, in any case, variability implies scatter and uncertainty. 



6 STATISTICAL METHODS 



There are two broad Galactic trends that relate to stellar ages: the age-metallicity 
relation (AMR) and the increase in disk scale height with age as stars and clusters 
encounter massive objects in the Galaxy. The AMR will not be treated here 
because its existence is in question but also because it is not very useful. For 
example, among stars fo rmed in the current epoch there is at least a factor of 
two range in metallicity. Wheeler. Sneden &: Truran (19891 ) have suggested that 
[0/H] may be a better chronometer than [ Fe/H] and that may be useful for 
populations but not for individual stars. See iFeltzing &: Bensbv (20091 ) for more 
on the AMR. 

The other phenomenon - known as disk heating - is gene rally studied as a 
means to understanding some Galactic dynamical processes (IWielen 19771 ). It 
can also be inverted to determine a mean age for a population if one has full, 
three-dimensional kinematic data. First, disk heating is an incompletely under- 
stood process and it is not clear what objects - giant molecular clouds, perhaps, 
or massive black holes - produce the observed effect. In any case, what is ob- 
served is a monotonic increase in the scale height of populations in the Galaxy 
perpendicular to the Galactic plane as one looks at progressively older ensembles 
of stars. Second, one needs full, three-dimensional kinematic information on an 
unbiased sample, preferably volume- limited. That means having good parallaxes 
and radial velocities in addition to proper motions. 

The classic reference on the subject is IWielen (mi\ ). who used Cepheids, 
nearby main sequence stars, and K and M dwarfs. The latter sample was di- 
vided by HK emission class, a rough age indicator available at the time. The 
main sequence stars were assigned ages to be half the MS lifetime for the mean 
of the sample. (This mean s that he did not correct the ages for the phenomenon 
that he was reporting on.) IWielen (19771 ) used the KM dwarfs because there were 
sufficient numbers of them close enough to the Sun to have good parallaxes at 
the time (before Hipparcos) . 

As an example, Soderblom (199(il ) used kinematic data to show that the BY 
Draconis binaries have a mean age of 1-2 Gyr and that the lack of old BY Dra 
systems was probably due to coalescence of the binary systems as they lost AM. 
Reiners Basri (20091 ) applied the method to a volume-limited sample of late- 



M dwarfs to derive a kinematic age of 3 Gyr. This is the same as the median 
age of G dwarfs in the solar neighborhood ( Soderblom. Duncan &: Johnson 199ll ) 
and results from the samples telling the same story: both G- and late-M dwarfs 
represent the full age range of the Galactic disk (~ 10 Gyr), and, given disk 
heating, the removal of the oldest stars from the solar neighborhood brings down 
the median age. 

The paper by Reiners &: Basri (20091 ) shows well how to apply the kinematic 
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method to a sample. 

7 THE LIFE CYCLES OF AGES 

Each phase in a star's Ufe offers specific problems and opportunities for estimating 
ages. 

7.1 Newly- Formed Stars 

The time-scales involved in the earliest phases of a star's life are 0.1 to 1 Myr 
or even less and so are very challenging to address. To a significant extent our 
assignment of the sequence of PMS phases is based on our preconceptions of how 
these systems ought to evolve, not actual measurements. The co-existence of 
several of these phases withi n a given star-form ing region is evidence that our 
understanding is incomplete. Evans et al. (20091 ) present a comprehensive look at 



the formation phases in nearby star-for ming regions and estimate relative time- 
scales and durations, and IWvatt f2008l ^ reviews the evolution of debris disks. 



Also pertinent for this age range is ai i alternate model for early accr etion that 
may affect ages at the ~ 0.1 Myr level (IWuchterl k Tscharnuter 20031 ). 



7.2 The Pre-Main Sequence 

Many of the most critical questions we ask about the formation and evolution 
of stars and planetary systems have to do with their lives before they reach 
the ZAMS, yet this phase presents severe challenges for determining ages. The 
most-used method is isochrone placement (Sec. 14. 2p . which is itself limited in 
every single aspect of its a pplication, from kn owledge of T^s and Lboi to the 
inadequacies of the models. iHillenbrand provides a review of age-datmg 



techniques for PMS star s using isochrones a nd other methods. Some other ap- 



proaches are provided by lMayne et al. (20071 ). who use empirical PMS isochrones 



to place clusters in an age order based on spline fits in magnitude versus color 
coordinates. This is a sensible means of applying the most basic information 
available, although the age ordering could be affected by differing age spreads or 
differin g binary frequencies or systematically erroneous distances, among other 
effects. iMayne et al. also try to place a "radiative-convective gap" in the 



CMD; this feature makes sense in terms of the physic s of models but in prac - 
tice appears prob lematic to locate precisely. Finally, Naylor &: Jeffries (20061 ) 



and iNaylor f;2009l ) describe a method for fitting the MSTO that takes account of 



evolution among massive stars. Their results yield ages that are 1.5 to 2 times 
those derived from PMS isochrones for very young clusters (~ 5 Myr old). This 
method is applied objectively but has the basic weakness of fitting the MSTO 
in a young cluster: a few of the most massive stars dominate the solution and 
their properties have the same vulnerabilities as other young stars (extinction, 
duplicity, etc.). 

Although it can only be applied to a group of stars, the detection of the lithium 
depletion boundary (LDB) is an especially promising technique for PMS ensem- 
bles in theoretical terms because the models used are straightforward. In practice 
the LDB method is of very limited use because it depends on obtaining good 
spectra of extremely faint objects, and that has been possible so far for only five 
nearby clusters (see Sec. 19. ip . 
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Asteroseismology may provide some help in dating PMS stars but will be lim- 
ited because our models of PMS stars are inadequate and because detection of 
oscillations in the presence of photometric variability is problematic. At present 
the relative faintness of most PMS stars puts them beyond the reach of the 
available seismic methods. However, more massive PM S stars lie in or near 
the instabilit y strip and can show detectable oscillations ( Guenther et al. 20091 . 



Zwintz 20081 ): these have been used more to test models of these stars than to 



establ ish ages. BD s in very young associations should also lie near the instability 
strip (ICodv 2009t ) and may reveal oscillations as well. 



The empirical methods should in principle work well for PMS stars because 
the changes in rotation, activity, and lithium abundances that occur in the PMS 
phase are large and rapid. But, at the same time, coeval populations of PMS 
and ZAMS stars in well-studied clusters show spreads in these quantities that 
exceed the slopes of the relations and there are saturation effects that lead to 
ceilings in these quantities for many stars. In other words, it is possible to apply 
the empirical methods if one has a group of 10-20 stars or more so that a mean 
relation can be determined, but the methods cannot be used reliably on individual 
stars. 

The movement of a PMS star in an HRD slows as it nears the MS, being roughly 
linear in terms of logr. Because of this there should be many more almost-on- 
the-ZAMS stars than there are T Tauri stars; these are the "post-T Tauris." The 
practical difficulty is that the post-T Tauris have had enough time to drift out 
of their formation regions into the field and thus are not so easily found. The 
major advance in this area in recent years has been the recognition that there 
are groups of PMS stars in the solar neigh borhood that can be iden tified on 
the basis of their activity and space motions ( Zuckerman &: Song 20041 ). Just as 



for the star-forming regions, these post-T Tauri ensembles often contain objects 
in a range of physical states, from classical T Tauri stars with active accretion, 
to "naked" T Tauris and post-T Tauris, emphasizing our poor overall understand- 
ing of these phases of stellar life. The TW Hydrae Association is a particular 
case in poi nt. The ages of these ensembles have bee n estimated from dyn amical 
expansion ( Mamajek 20051 ) and lithium abundances ( Mentuch et al. 2008 ). In all 



cases it is the average age of the ensemble that is determined. 
7.3 Zero-Age Main Sequence and Near-ZAMS 

As is well known, stars change very little in their structure for the first ~ 1/3 of 
their MS lifetimes; this is why they form a well-defined main sequence, after all. 
For massive stars, their MS lifetimes are brief and, in part because of that, they 
are almost always found in star- forming regions and associations. In that case the 
ages can be estimated either from the MSTO of more massive stars, or from where 
less- massive stars are just reaching the ZAMS. The practical problems in these 
cases are those of open clusters (see Sec. [9|). Isolated massive stars are rare but 
do exist; estimating their ages is difficult because the luminosity is uncertain, but 
it is sometimes possible when a spectroscopic gravity indicates the evolutionary 
state. 

Intermediate-mass stars (late B through mid-F) can be difficult to age-date if 
they are isolated and removed from a context provided by a cohort. However, A 
and F stars in particula r have had t emperatures and gravities determined from 
Stromgren photometry ( Olsen 19941 ). These can be calibrated against nearby 
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open clusters which span an age range comparable to the MS lifetimes of A 
and F stars, but the uncertainties are again significant at least in relative terms. 
Recent advances in stellar interferometry, however, have shown that some A stars 
are measurably non-spherical, meaning that traditional one-dimensional models 
cannot predict their surface properties and that we may have little idea of their 
evolutionary state. 

Solar-type stars on and near the ZAMS are discussed elsewhere in this review. 
Lower mass stars (late K and M), however, remain challenging. The empirical 
methods that work for G dwarfs are less helpful at lower masses. Getting an age 
from a rotational period is still effective, but lithium is depleted far too quickly to 
be of use in most cases. Ac tivity in the low est-mass stars is notoriously variable. 



even for old stars, although I West et al. (2008 ) have studied activity as a function 



of age in M dwarfs by using a large sample and statistically deriving ages of 
subsets from their kinematics. 

7.4 Subgiants 

Subgiants are the "sweet spot" in the HRD for determining ages from isochrone 
placement. These stars are clearly well-removed from the main sequence and 
are changing quickly in luminosity and T^s- These ages are obviously model- 
dependent, and metallicities are needed to know which set of models to which 
to compare, and the fundamental limitations are in the models themselves and 
our knowledge of the physics of such objects. Intermediate-mass stars and lower 
masses develop or deepen CZs, and this can alter surfaces abundances of some 
elements. Stromgr en photome try can be used to determine T^g and log (7 for A, F, 
and early-G stars ( Olsen 19941 ) to place them in HRDs, but relative uncertainties 
are significant. 

7.5 Very-Low-Mass Stars 

Solar-type stars age slowly but very-low-mass (VLM) stars evolve hardly at all 
once they reach the ZAMS. However, like all stars, prior to the ZAMS they change 



significantly in Lboi (jBurrows et al. 200ll ) as do brown dwarfs (BDs) . All of these 



objects that are near the bottom of the stellar mass function are of interest in 
part inherently, in part for their role in the Galactic mass budget, and in part 
because they form the transition from bona fide stars to planets. The obvious 
way to limit the age of a VLM object is to associate it with another star or 
group for which the age is better constrained. Lacking that, these objects are 

c hallenging. 

Jameson et al. f;2008l ) describe a method to determine ages of L dwarfs from 



near-infrared photometry and distances when those objects are younger than the 
Hyades (i.e., < 0.7 Gyr). The method uses Mk, the absolute iC-band magnitude, 
and is c alibrated against nearby c lusters. Their estimated uncertainty in logr is 
0.2 dex. iBurgasser fc Blake (20091 ) were able to place a lower limit on the age of 



a binary system composed of an M dwarf with a T dwarf, again because having 
a binary constrains models significantly. Howe ver, association of a V LM object 
with another group is not always the solution: Bonnefov et al. (20091 ) show that 



a VLM binary thought to be a member of the ~ 10 Myr-old TW Hya association 
is itself about 30 Myr old and so unlikely to be a bona fide member. 
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Burgasser f 20091'! lists five age es timation techniques for VLM objects and com- 



pares them. iDupuv. Liu Sz Ireland (2009 ) illustrat es several case studies applying 



these methods to brown dwarfs. West et al. (20061 ) studied Ha activity in a large 



sample of Sloan Digital Sky Survey M7 dwarfs and were not able to delineate an 
activity-age relation but found those low-mass stars decline in Hct rapidly at an 
age of 67 Gyr. 

7.6 Very Old Stars and Population II 

Most of the Galaxy's oldest stars are too distant to currently have good paral- 
laxes which would allow the luminosity to be calculated to place the star in an 
HRD. Instead it is possible to use i iucleocosmochronome try (Sec. 13. ip . Another 
possibility is using Be abundances ( Smiljanic et al. 20091 ). similar to using Li for 



Pop I stars. 

8 CALIBRATING AGE INDICATORS 

The most fundamental problems in applying age estimators to stars arise from 
the difficulty and uncertainty of calibration. Two methods will be discussed here: 
using clusters to define a property at specified ages, or taking advantage of having 
a volume-limited sample. Both methods require assumptions to be made, and 
those assumptions can be difficult to test. 

8.1 The Dissolute Lives of Open Clusters: The Sun as Anchor 

The obvious way to calibrate an age relation is to observe stars in clusters. This 
still leaves the problem of establishing the ages of clusters, but clusters have 
many stars and so it is usually possible to determine mean properties with 
good precision. The underlying astrophysical limitation is that basic Galao 



tic forces tear PCs apart o n a time-scale of ~ 200 Myr (jjanes &: Phelps 1994 : 



Janes. Tilley Sz Lynga 19881 ) and more weakly-bound clusters can simply fall 



apart from their own internal motions more quickly than that. Old open clus- 
ters are inherently rare and not as near as younger examples. Making matters 
worse, the age-related quantities to be measured (rotation, activity, and lithium 
in particular) in the older stars get progressively weaker, requiring ever higher 
signal-to- noise spectra and high resolution for ever fainter stars. 

This process of cluster dissolution, of course, is where field stars come from, 
and so the stars we most want to place ages on are least represented among 
potential calibrators. There is also a potential bias in calibrating against the 
few old clusters we can observe. An old cluster such as M67 was probably born 
unusually rich and dense to have survived in the Galaxy for so long. As such, 
its properties may not be fully representative of field stars of similar age, which 
came from more ordinary backgrounds. In particular, the partitioning of angular 
momentum within a dense cluster could be very different from that for a loose 
association and that could have ramifications even late in life for stellar spindown. 

On the other hand, among the older stars we have the Sun, for which many 
properties are precisely defined, and so the Sun is generally used to anchor age re- 

uestion of whether 



lations using its behavior at 4.5 Gyr. This naturally raises the q 

the Sun is itself representative of a 4.5 Gyr-old 1.0 Mq star. Gustafsson (20081 ) 



(and others) have argued that the Sun is peculiar or unusual in one respect or 
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another compared to other G dwarfs. But given the many well-defined properties 
that characterize the Sun, it would be odd if several did not deviate significantly 
from statistical means, just as each of us is both typical and unusual as a human. 
The null hypothesis should be that the Sun, on the whole, is a typical star for 
its mass, composition, and age, and that hypothesis has not been convincingly 
disproved. 



8.2 Clusters as Benchmarks 

Age determinations for OCs will be discussed in Section [9l They clearly provide 
the primary benchmarks we use to study age-related properties, not just to define 
a mean for a given age but also to examine how well a coeval population behaves. 
Also, as noted in Sec. 14.2.21 the available evidence suggests that any age spreads 
within clusters must be small, ~ 1 Myr, and so the dramatic scatter we see in 
various properties among stars in ZAMS clusters must be explained in terms of 
the phenomenon itself. 



8.3 Volume-Limited Samples 

Volume-limited, or at least unbiased samples provide a special opportunity be- 
cause all ages must be fairly represented. That being the case, other criteria can 
be used t o study ages and age relations, although assumptions must be made. For 



example, ISoderblom. Duncan &: Johnson (199ll ) studied a volume- limited sample 



of activity in G dwarfs and were able to derive an age-activity relation by using 
the Sun to anchor 4.5 Gyr and assuming that activity declines monotonically with 
age (or at least that the decline with age is dominant compared to variability), 
and then also applying a correction for "di sk heating" to account for the lack 
of very old stars in the solar neighborhood. iReiners &: Basri Hooi) have used a 



volume-limited sample of late-M dwarfs to determine a mean age by calculating 
their Galactic kinematics and comparing to observations of disk heating (see Sec. 

ED. 



8.4 The Special Role of Binaries 

Binary systems can be especially helpful in calibrating and checking age relations. 
First, binary components provide a consistency check in that both stars should lie 
on a single isochrone (to within the errors) for a given age indic ator if that indica- 
tor can be measured in both stars. For example, Barnes (20071 ) used rotation ob- 



servations of well-separated visual binaries to co nfirm his gyrochronology calibra. - 



tion. HK activity fares less well in the study of Mamajek Sz Hillenbrand (20081 ) 
when components in binaries are compared, but most systems meet this expec- 
tation, even if there are conspicuous exceptions. 

Second, binaries provide "mini-clusters." A system with, say, a solar-type sec- 
ondary and an intermediate-mass primary has its a ge limited by the HRD location 



of the more massive star. This method was used by Soderblom. Duncan &: Johnson (199ll )| 



in calibrating HK emission versus age. The obvious limitation is that a binary 
provides just a single datum, nowhere near the ensemble average one can achieve 
for a cluster. At the same time, a system such as a Centauri provides not just 
well-studied stars but well-determined masses, and that is rare. 
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Binaries can also be helpful on setting l imits to an age when one star i s sig- 
nificantly more massive than the other. iGahm. Ahlin &: Lindroos and 
Lindroos (19831 ) . for example, sought F, G, and K dwarf companions to O and 



B stars as a means to identifying potential post-T Tauri stars, solar-mass stars 
that are ~ 10 Myr and difficult to find in the field. However, in many cases 
these necessarily wide pairs, as judged by som e of the age criteria discussed here , 
were optical doubles, not physical b inaries (iMartm. Magazzri fc Rebolo 199^ : 
Pallavicini. Pasquini fc Randich 19921). although different studies do n ot agree 



on which systems may be physical ( Gerbaldi. Faraggiana &: Balin 200ll ). 

A va riation on binaries as mini-clusters occurs when one of the stars is a white 
dwarf. Catalan et al. (200^ ) have used WDs in wide binary systems as a check on 
ages of low-mass stars. Wide binaries (with separations of ~ 1000 AU or more) 
are used as te st particles for age-related dyn amical processes, both in star-forming 
regions (e.g., iKraus k. Hillenbrand f 2009^^1 1 and amon g the older populations of 
the disk (e.g., Makarov. Zacharias &: Hennessv (20081 )). For those systems one 
can help weed out non-physical pairs with independent age-related criteria. 

Other special uses of binaries are discussed in the next section on OC ages: 
eclipsing binaries as mass calibrators and tidal effects as an age-ordering method. 

In some instances binary systems may be misleading. As just noted, one must 
first apply astrometric and radial velocity criteria to ensure the two stars form 

a physical pair. But even when they do the masses we see now are not necessarily 

the same as when the system formed. For example. iRappaport. Podsiadlowski &: Horev f2009l l 
show that the primary in the Regulus system has a significantly different mass 
now than it did in the past; Regulus is one of the Lindroos systems just noted 
and this fact, of course, substantially alters conclusio ns one might draw about 
the companion K dwarf. In another potential example, IStassun et al. f2008l ^ note 
that the two components of the PMS system Parenago 1802, in the Orion Nebula 
Cluster, have the same mass yet substantially different temperatures and lumi- 
nosities. They suggest that the two stars may have formed at different times 
and that is possible, particularly if companions may have been exchanged as two 
binaries encountered each other. However, this system has an orbital period just 
under 5 days which would make it resistant to dynamical effects but, at the same 
time, the two components of this system in their earliest stages were significantly 
larger than they are now and their propinquity may have led to effects that pro- 
duced the curren t differences in Tpff and Lbni • Simila, rlv, the brown dwarf eclipsing 
binary found by IStassun. Mathieu &: Valenti (20071 ) may be in a different state 
now than when it formed; its orbital period is just under 10 days. Things are not 
always what they seem. 



9 OPEN CLUSTER AGES 



OCs provide the benchmarks to calibrate virtually every stellar age except that 
of the Sun. The essential tool for studying clusters is the HRD - log Lboi ver- 
sus logTgff - and its observational counterpart, the CMD. HRDs are powerful 
tests of stellar physics because the models have to work consistently over a very 
broad range of mass and evolutionary states, where many different physical pro- 
cesses are taking place. At the same time, that means that derived ages to 
some extent themselves depend on the cluster age. For instance, very young 
clusters usually have ages determined from their massive stars at the MSTO, 
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whereas intermediate-age clusters have ages determined from intermediate-mass 
stars where the interior physics is different. Old clusters' ages often depend on 
the locations of evolved stars in the HRD. 

For the present purposes, the relevant OCs are the nearby well-studied exam- 
ples (within about 1 kpc). The more distant clusters are important for under- 
standing our Galaxy, particularly because old OCs are inherently rare, but it 
is primarily the nearby OCs that will be discussed here because their low-mass 
members are bright enough to be accessible to spectroscopy and with that a 
fuller understanding of their properties. In most cases there is ample observa- 
tional information to assess the membership of individual objects well, reducing 
systematic effects from field star contamination, and the metallicity of the cluster 
is known from spectrum analyses. These are some of the problems and limitations 
that are considered in analyzing the ages of the nearby OCs: 

1. Reddening and extinction are always uncertain, but for the nearby OCs 
those quantities are small and so the errors tend to matter less. However, 
PMS clusters and associations are generally still located in and amid the 
material from which they formed and their reddenings can be high and 
significantly different from star to star. The reddening law to apply also 
may vary from place to place and is a source of systematic uncertainty. 

2. Both Lboi and Tefj are theoretical quantities that are useful for describing 
models, but are simplifications when considering real stars. In particular, 
Teff is a parameter that describes the energy flux of a star through its sur- 
face, based on calculations of the energy generation occurring in the stellar 
interior, but Tefj as a single number does not describe the way that flux 
appears in observations. Stellar surfaces are complex and the formation of 
spectrum lines and spectral energy distributions are also complicated when 
considered in detail. A star like the Sun, for example, would need three com- 
ponents to describe its surface even in a fairly simplistic way: the upwelling 
hot material, the cool material moving downward, and the spots, each with 
an areal coverage and a temperature or temperature gradient. Chromo- 
spheric activity also changes the temperature structure in the outer parts 
of a star's atmosphere, altering the line formation conditions. It is this in- 
herent complexity that likely accounts for how different photometric indices 
predict different temperatures, none of which agree with what comes from 
spectrum analysis. Comparisons of different color-temperature relations 
(CTRs) with each other and with Tp ff values determined from spectrum 
analysis show differences of 80-100 K ( Ramirez fc Melendez 200*51 ) and 100 



K is generall y taken to represeii t the true uncertainty of stellar Tgfj values. 
The study of lRibas et al. (20091 ) is illustrative in that they derive T^s val- 
ues for a single, well-observed star using a variety of methods (photometry, 
ionization balance, excitation balance) without flnding a single, consistent 
value. A part of the problem is that different color indices are sensitive to 
different aspects of a star's atmosphere. In particular, (B — V ) appears to be 
sensitive to spots on young G-K dwarfs (jStauffer et al 2003ll and the (V-I ) 



index appears to produce more reliable T^fi values ( Stauffer et al. (20031 ): 
Lvra et al. (200dl )). 



3. OC ages derived from different methods with different underlying physical 
processes can yield ages that different systematically to a significant degree. 
A notable example are ages from the lithium depletion boundary (LDB) 
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in very young clusters and groups (see below). Five LDB ages are now 
available and all are systematically about 50% larger than ages determined 
from the MSTO. The age scale from the LDB involves fundamentally simple 
physics, whereas MSTO ages rely on knowledge of the interiors of massive 
stars where processes may be taking place that we do not fully appreciate. 
Because of this, the LDB age scale may be the more reliable, but there are 
too few LDB ages to recalibrate the overall OC age scale (see Sec. 19. ip . 

4. The transformation of magnitudes to logLboi requires knowledge of the 
cluster's distance, the extinction, if any, and bolometric corrections. Only 
a very few OCs have trigonometrically determined distances, and even some 
of those remain contentious. Indeed, the process of fitting isochrones to OC 
CMDs generally yields a best-estimate distance as part of the exercise. 

5. Because the initial mass function for clusters is strongly biased against 
high-mass stars, there are generally very few stars in young clusters to 

define the MSTO and this is exacerbated by the p resen ce of binaries. 

■Torgensen fc Lindegren (20051 ) . iNavlor fc Jeffries (200(^) and lMavne fc Navlor r2008l ) 



have attempted to reduce systematic errors in MSTO fitting by applying 
bayesian techniques. Similar efforts for many more of the nearby clusters 
could at least lead to a more consistent age ordering. 

Mermilliod (2006 ) reviewed age estimation methods for OCs and included some 



of th e techniques that use HRD m orphology that are mainly a pplied to old clus- 



ters. iGallart. Zoccali fc Aparicio (2005 ) and .Mevnet et al. (200a ) provide a thor- 



ough discussion of the problems faced in determining cluster ages from isochrones, 
particularly a close examination of stellar models. That review is oriented toward 
disentangling the star formation histories of nearby galaxies from their composite 
CMDs, but the limitations and problems that arise in applying the models are 
discussed extensively. For example, an age derived from a MSTO can change by 
up to 50% depending on the specific models used, but that difference is greatest 
in certain age ranges and very low in others. This is a reflection of the different 
non-standard physical processes that may be included by different groups because 
those processes take place within specific mass ranges, which is to say age ranges, 
given the obvious connection between the masses of s tars at the MSTO and the 
clus ter age. For so me particular well-studied cases, Pinsonneault et al. (20041 ) 



and An et al. have been able to constr uct empirical isochrones for ZAMS 



and post-ZAMS OCs, and lMayne et al. r2007l ) have done so for PMS stars. This 



was done to lessen problems with the CTRs and led to improved distance esti- 
mates but did not directly affect derived ages as much. 

The physical processes that are not well understood include convective core 
overshoot, rotationally-induced mixing, internal gravity waves, and diffusion. 
Magnetic fields may also play a role. This review cannot treat this subject 
in any detail, but what is important here is that several of these processes - 
convective core overshoot and rotationally-induced mixing in particular - replen- 
ish fuel in the nuclear-burning regions of stars a nd so can significantly extend 
the MS lifetimes of mas sive stars by up to ~ 50% ( Rosvick &: VandenBerg 19981 . 



Woo fc Demaraue 200lh . This leads to fundamental and significant uncertainty 



in the overall age scale of clusters, especially young clusters, and of stars in gen- 
eral. 

Very young clusters present their own challenges. For any method, the un- 
certainties in both Lbol find T(,s are significantly larger for PMS stars and this 
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has many causes, some observational and some inherent t o the stars themselves; 
see Sec. 14. 2i Along the upper MS, iMavne et al. f2007l l discuss a "radiative- 



convective" gap where stars develop their radiative core as they first reach the 
ZAMS; this gap is very close to where the PMS line for low- mass stars meets the 
ZAMS. The R-C gap moves in the CMD with age and allowed iMayne et al. 



to set an age ordering for very young clusters that is different from the age order 
determined from the MSTO. However, the R-C gap is a lo w-contrast feature in the 
CMD that can be difficult to perceive in sparse clusters. iMavne fc Navlor ^20081 ) 



applied objective criteria ( Navlor &: Jeffries 20061 ) to MS fitting of very young 



clusters to avoid the biases that come with fitting by eye. A major goal of these 
studies is to determine the extent of age spreads within star-forming regions, and 
that will be discussed in the next section as well. 

PMS clusters can sometimes have conflicting age indicators, as when massive 
stars a re leaving the MS at the MSTO and low-mass stars are just reaching the 
ZAMS. Lvra et al. (20061 ) have found that such discrepancies depend on the color 



index used, as noted above, and that the MSTO and MS turn-on ages are in good 
agreement when colors such as {V — I) are used that avoid the blue part of the 
spectrum. 

Old clusters also pose difficulties. Even for well-studied cases such as M67 and 
NGC 188 there can be field star contamination among the faint stars of the lower 
main sequence and the degree of convective core overshooting remains uncertain 



(jSalaris. Weiss k Percival 20041 : IVandenBerg k Stetson 200' 



9.1 The Lithium Depletion Boundary 

The potential significance of these alterations of the OC age scale is illustrated by 
the ages derived from the location of the lithium depletion boundary (LDB) seen 

in some young clusters. The LDB was first proposed as a means of discriminating 

bona fide brown dwarfs (BDs) from very-low-mass stars by lRebolo. Martin fc Magazzii 



They pointed out that objects below about 0.06 Mq will never get hot enough in 
their cores to destroy their Li, which happens at temperatures above about 2.5 
MK. The temperature and luminosity of that dividing line changes with age as 
the 0.06 Mq stars cool and as less massive objects acquire the core temperature 

n eeded to astrate Li. 

Basri. Marcy &: Graham were the first to detect the LDB in a clus- 



ter. The y detected Li in a single very faint Pleiades object and applied the 
models of Nelson. Rappaport &: Chiang (19931 ) to derive a cluster age of 110-125 



Myr, substantially greater than the 7 Myr that results from fitting the MSTO 
Stauffer. Schultz Sz Kirk patrick 1 9981 ). Subsequent work has refined the models 



Burke. Pinsonneault Sz Sills 2004 ) and has added new observa tions that bring 



the numb er of detections of th e LDB to five clusters: the Pleiades (IStauffer. Schultz fc Kirkpatrick 199 



and IC 4665 (|Manzi et al. 20081 ). 



a Persei (IStauffer et al. 1999ll: N GC 2547 and IC 2391 (jjeffries k Oliveira 20051 ) 



For the youngest clusters the absolute differences in age between MSTO and 
LDB determinations are small, yet the LDB ages in all cases are systematically 
about 50% greater than MSTO ages. The observations needed to detect the LDB 
in a cluster require good spectra of extremely faint objects, limiting the test to 
the nearer OCs. Yet the LDB ages pose a real problem, and that is because the 
physics in the models that predict the position of the LDB with age depends only 
on the core temperature reached and that is calculated straightforwardly. The 
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surfaces of BDs may be more complex than we can understand at this time, but 
that doesn't influence what occurs in the object's core. For these reasons the age 
scale implied by the LDB observations may be more reliable than that from the 
MSTO, but there remain too few clusters with LDB ages to readjust the scale of 
MSTO ages, and that situation is not likely to change for some time. 

It is also possible, however, that our estimates of the radii of very young 
stars are systematically wrong, and this would affect the interpretation of LDB 
ages. There is accumulating evidence that PMS and ZAMS stars are larger tha n 
models predict fe crger et al. (2006); Lopez-Morales (2007); Rib as et al. (20081 ): 
Torres. Andersen k Gimenez (20091 ): I Jackson. Jeffries k Maxted (20091 )) and thisl 



makes them appear cooler; adjusti ng for this effect ma y help to bring LDB and 
MSTO ages of OCs into agreement ( Yee k Jensen 2o"lol). Other evidence suggests 



the MSTO ages are sound. For instance. iMakarov (20061 ) was able to precisely fit 



a model to the star a Per itself that is consistent with a good fit to the MSTO in 
the « Persei cluster fo r an age of 52 Myr, as compared to an LDB age of 90 it 10 



Myr (jStauffer et al. 1999i ). 



9.2 Other Indicators of Cluster Ages 

9.2.1 TIDAL EFFECTS ON BINARY ORBITS Tidal interactions be- 
tween the stars in close binary systems influence the evolution of the stars them- 
selves if they are close enough but more generally lead to changes in the orbit, in- 
cluding synchronization of the stellar rotation with the orbit and circularization of 
the orbits. This is especially effective in low-mass stars because they can dissipate 
angular momentum in their winds and thus remove AM from their orbits, leading 
to coalescence. The theory of these interactions is not important here; what mat- 
ters is that we expect to see a progression with cluster age of the period at which 
synchro nization and circularization are seen and that is largely borne ou t observa- 
tionally (jMeibom k Mathieu 20051 . iMeibom. Mathieu k Stassun 20061 ). although 



the Hyades and Praesepe appear to conflict with the trend. 

Whether or not the physics of tidal dissipation is fully understood, the process 
can only act in one direction and so observations of these tidal effects should make 
it possible to place OCs in a reliable age ordering. In practice the observations 
required are intensive (precise photometry and spectroscopy over extended times 
to get rotation periods and orbital parameters) and limited to nearer clusters. 

9.2.2 WHITE DWARFS IN CLUSTERS The detection of the white 
dwarf cooling sequence in a cluster allows for a determination of the cluster's ag e. 
This subject has been we ll reviewed elsew here ( Koester (20021 ). [von Hippel (200 5) .[ 



Wineet k Kepler (20081 ). ISalaris (20091 )) and so only a few comments will be 



made here. 

First, the physics of WD cooling itself seems straightforward because no energy 
is being generated; there is just gradual cooling of the star with time and the 
physics of that process is fairly well understood. However, there are uncertainties 
that matter for the input physics, such as in the equation of state, opacities, etc. 
The detailed energy budget is strongly dependent on the chemical stratiflcation of 
both the degenerate core and the envelope. All of these effects (and more) mean 
that some claimed uncertainties for cluster ages from WD cooling sequences are 
probably optimistic, and there is an inherent error of 1 Gyr. 

Second, to derive the cluster age also requires determining the initial-flnal mass 
relation, which relates the WD mass to the mass of the progenitor star prior to 
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mass loss on the asymptotic giant branch. There are inherent uncertainties in this 
determination, such as the degree of mass loss on the asymptotic giant branch. 

Third, the above uncertainties arise from the models, but there are additiona l 
observational factors (including composition and temperature); see I S alar is 



Fourth, the usual technique is to detect the coo lest WDs in a cluste r but that 
can be difficult and can lead to selection effects. iJefferv et al. discuss a 



means of using brighter WDs in clusters to the same end. 

Finally, the cluster as seen today can be the result of dynamical evolution that 
influence s the conclusions drawn . An example is the WD age determined for the 
Hyades (IWeidemann et al. 1993). whi ch, at 300 Myr, is about half the MSTO 



age of 625 Myr (jPerrvman et al. 19981 ). This discrepancy has been argued to be 
due to the coolest WDs having been ejected from the cluster. 
9.2.3 ECLIPSING BINARIES E clipsing binaries (EB s) play a special role 



in studying stellar astrophysics and OCs ( Stassun et al. 20^ ). If the components 



are well separated so that mass exchange or other physical interactions are un- 
likely to have taken place, then finding EEs in OCs and measuring their orbits 
allows one to place a known mass benchmark on the models. In addition, EE or- 
bits yield radii, another fundamental prop erty of the models. As a resul t EEs can 
provide critical tests of model isochrones ( Southworth Sz Clausen 20061 ). particu- 



larly in the fortuitous case of an EE being found near the MSTO. Also, EEs al- 
low accurate distan ces to be determined, even for nearby galaxies ( Clausen 20041 . 



Guinan et al. 19981). and those distanc es are reddening- and extinction-independent .| 



For instance, iGrundahl et al. (200811 studied an EE in the metal-rich old clus- 



ter NGC 6791, and in NGC 188 iMeibom et al. (20091 ) have analyzed an EE near 



the MSTO. In both cases it was possible to find a precise age (to about 2-3%) for 
the EE by comparing to theoretical isochrones in a mass-radius diagram, inde- 
pendent of distance, reddening, and a color-temperature relation. This enabled 
them to critically test conventional isochrones and their sensitivity to metallic- 
ity, distance, and reddening. The ages that result from analyzing EEs are still 
model-depende nt, but offer an addi tional means of testing stellar physics. In a 
similar fashion. IClausen et al. Hooi) were able to use models to get an age for the 



field EE V636 Centauri, a pair of solar-type stars. This was especially important 
because the radii for this system were at odds with conventional models, which, 
along with other evidence, has suggested that young, active solar-type stars have 
larger radii than has generally been assumed (see comments in Sec. 19. ip . 

9.3 OC Age Precision and Accuracy 

For reasons noted above there remains significant uncertainty in the accuracies of 
OC ages. Even ignoring discrepancies between different m ethods, the standar d 
technique - isochrone fitting - has a precision of ^ 10% ( Mevnet et al. 20091 ). 



As an example, the Hyades presents what has to be the most favorable case 
possible, with zero reddening, a well-determined metallicity, only a very small 
uncertainty in distance, a fairly rich population, and unambiguous membership 
and multiplicity. Even with these significant adva ntages the best age is 625 ifc 50 
Myr ( Perrvman et al. 19981 ). As another example, J0rgensen k Lindegren (200^ ) 



applied a bayesian fitting method to the MSTO in two clusters: IC 4651 (deriving 
r = 1.56 ± 0.03 Gyr) and M67 (r = 4.05 ± 0.05 Gyr). These are remarkably 
small errors (~ 1 — 2%) but do not take into account the full range of possible 
uncertainties. The uncertainties usually quoted for OC ages are ~ 10% and that 
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probably represents a good estimate of just the fitting errors. 



10 A COMPARISON OF TECHNIQUES AND A SUMMARY 



How well do results from different methods compare? First, we note that some 
methods are inherently limited to being applied to ensembles of stars, while 
others are for individual stars. For individual stars, the available methods include 
isochrone placement; nucleocosmochronometry; asteroseismology. The methods 
limited to ensembles include: 

• The lithium depletion boundary (LDB) method for young clusters, because 
it depends on finding the point among a group of low-mass objects where 
lithium reappears. 

• Kinematic expansion, because it involves determining when in the past a 
group of stars occupied the least space. 

• The empirical methods - rotation, activity, and lithium - when applied 
to young stars (r < 0.5 Gyr) because the inherent scatter among coeval 
populations is large and exceeds the change in the quantity over significant 
times. 

• Isochrone fitting to OCs, for obvious reasons. 

And the empirical methods (declines in rotation, activity, and Li) can be used 
for either individual stars that are > 0.5 Gyr or ensembles that are younger. 

There are not many instances of stars with ages determined in multiple ways 
so that one can compare results from the various methods. Also, methods such 
as using rotation and activity are not really independent of one another and show 
an inherent physical correlation. However, Figure [8] compares ages of G dwarfs 
derived in two very different ways: from isochrone placement (Sec. 14. ip and from 
rotation rates (Sec. 15. ip . Th e lack of a perceptible correlation between the two 



of nearby main sequence stars using available information (isochrone placement, 
rotation, HK emission, space motions, and metallicities) and similarly saw poor 
agreement among the indicators. 



If this review has seemed downbeat on the prospects of determining ages of 
individual stars that is because there is a lot to be negative about and I have 
wanted to ensure that readers appreciate how difficult the task is. Few topical 
symposia have been held on the subject of ages; the first was "I'Age des Etoiles," 
held in Paris in 1972. A look throu gh the proceedings of that lAU colloquium 
(ICavrel de Strobel fc Delolace WT^ ) shows that many of the subjects discussed 
there have advanced rather little. We now have much more and much better 
observations, but many of the same problems remain. 

Also, one desires to determine an age with a precision comparable to other 
stellar quantities and there are inherent factors that prevent that. Nevertheless, 
age-dating is not an impossible task, just a slippery and uncertain one, and it is 
possible to establish at least limits on how old most objects are. Despite what I 
have said, I am very optimistic about the future. How can we improve on matters, 
and what does the future hold? What can be done to improve the situation? 



Fig. [8l HK ages ver- 
sus ages from isochrone 
placement. 




derived ages for a number 



11 THE FUTURE OF AGE 
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First, the near future is very promising, in large part because of the CoRoT and 
Kepler missions. In particular, Kepler is obtaining ultra-high-precision (1:10^) 
photometry for ^ 100, 000 solar-type stars in order to detect the transits of 
Earth-sized planets. In addition to what it does for planet-hunting, Kepler will 
leave a rich legacy for stellar astrophysics, for never before have we been able to 
observe stars at such a precise level (excepting the Sun itself, of course, and the 
advent of similar quality photometry for the Sun is very recent). The Kepler data 
will be like "full body scans" for solar-type stars and should enable the detection 
of not just rotation in older stars but also differential rotation as the apparent 
Prot changes. 

Of greatest relevance here is the fact that Kepler obtains its photometry at 
one-minute intervals for several thousand selected, brighter stars in its field (the 
remainder are sampled in 30-minute integrations). High photometric precision 
combined with long on-target time sequences will enable the detection of solar- 
like oscillation modes in these stars. Those asteroseismological detections will 
enable reasonably precise ages to be assigned, and most of those stars will be 
older stars of the field. Those ages will only be as good as the models used 
to calculate stellar properties, but, at the same time, the physical constraints 
imposed by asteroseismology will stress those models and help us refine them. In 
addition, Kepler'' s high precision should make it possible to detect, for instance, 
those stars' rotation periods, even for the old, inactive stars. The apparent P^ot 
changes as spots migrate, as they do on the Sun, and so good photometry can 
also enable us to measure stellar differential rotation with precision, a critical 
parameter for dynamo models. With sufficient devoted time on selected targets 
(months), asteroseismology can also show the signature of the convection zone 
depth, the first time we will have such a key constraint on models for any star 
other than the Sun. Kepler is a solar physics mission at heart. These ages from 
asteroseismology for older stars will be a breakthrough in that they will enable 
much better testing of the empirical methods and should help us understand if 
ages from isochrone placement and the empirical methods are sound. 

Another key mission in the near future will be Gam, which will obtain highly 
precise [a ~ 100/iarcsec) parallaxes for millions of stars in our Galaxy. Gaia 
should remove almost all ambiguity about the true distances to any of the open 
clusters used as calibrators and will also help weed out non-members in those clus- 
ters. With results from Gaia, luminosities for the nearer stars will be virtually 
error-free except for uncertainty in bolometric corrections. Gaia will also pro- 
vide accurate distances to PMS objects in star- forming regions, again removing 
ambiguity and helping us to understand those stages of stellar evolution. 

The proposed variability surveys, such as the Large Synoptic Survey Telescope, 
will provide a bounty for studies of stellar variability. More important, many 
eclipsing binaries will be discovered, enough so that we can choose those best 
suited to test stellar models. In general these missions improve our knowledge 
and understanding of our Galaxy by probing the data domain along dimensions 
not yet explored. 

Our ability to study, understand, and model the physics of stars has improved 
enormously over my career. The new domains of astronomical observation that 
are becoming available will only challenge what we think we already know and 
will lead to substantial improvements. Each new domain we explore reveals the 
weaknesses of the models and forces improvements. 
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11.1 Predictions 

If I were any good at predicting you would have heard about it before now. 

1. Kepler will produce excellent asteroseismological data for a large number 
of solar- type stars and that, combined with additional observations and 
careful modeling, will yield better ages for older stars than we have ever 
had. Those ages will enable us to calibrate and understand empirical age 
indicators such as rotation and activity. Rotation will turn out to be as 
well-behaved as we have assumed. Activity will be useful statistically for 
stars more than ~ 1 Gyr old, and we will gain some knowledge of the 
inherent variation and scatter in activity at different ages. Lithium will 
remain confusing and resistant to modeling. Kepler will also for the first 
time turn asteroseismology from a boutique operation into a wholesale one, 
with good frequencies available for thousands of stars. 

2. Critical physical information, such as the depth of stellar convcctive en- 
velopes, will result from long on-target integrations with Kepler. When 
combined with concurrent observations of chromospheric activity we will 
gain insight into the factors that drive dynamos and long-term activity cy- 
cles. If we can understand the physics of the solar 11-year cycle we can 
start to understand stellar dynamos in general and we can turn rotation 
and activity from empirical age indicators into model-dependent ones. 

3. Suitable binaries will be found in young clusters enabling accurate masses to 
be measured for ZAMS solar-type stars. Those masses will be significantly 
different than the models have predicted for a given location in the HRD 
and will lead to important insights in stellar models. This may hold the key 
to understanding the scatter seen in rotation, activity, and lithium in young 
clusters, as well as the discrepancy between ages from the main sequence 
turn-off and the lithium depletion boundary. The large scatter seen for 
rotation, activity, and Li in PMS and ZAMS clusters will be explained as 
being due to the distribution of angular momentum, disks, and companions 
when stars are formed and the later consequences of those. 

4. The best-available data for the nearest OCs will be analyzed and fitted in a 
self-consistent manner to produce at least a more reliable age ordering for 
those clusters. This will lead to some insights on what factors produce dif- 
ferences between clusters that we do not now understand, such as conflicts 
between MSTO ages and an age determined from PMS isochrones. 

5. Knowledge of age and mass will enable us to create a sample of bona fide so- 
lar analogs and we will see that the Sun is a very typical star. Accurate ages 
will make it possible to find a solar-age star in a "Maunder minimum" state 
and that will help us understand better the connection between changes in 
luminosity and activity for stars like the Sun. 

6. An Earth-like planet will be found around a nearby solar-type star and 
that planet will show a biomarker. We will estimate the star's age from 
asteroseismology and it will be > 2 Gyr old. 
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13 Sidebar 1: Why Care About Stellar Ages? 

Some day soon we anticipate the detection of Earth-sized planets around stars. 
Before too long, we hope, we can expect that someone will identify a sign of life 
— a biomarker - on an Earth-like planet around another star. Given the need 
for high angular resolution to make that observation and the inherent paucity of 
photons to work with, the star in question will be one close to us, which is to say a 
field star that is unassociated with a cluster. When that claim is announced, the 
first question we will all ask is " How old is that star?" , so that we can assess the 
planet's evolution. Yet determining that age is likely to be difficult and imprecise, 
as problematic in its own way as the observation of the biomarker. 
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14 Sidebar 2: A Challenge for Stellar Astrophysics 

For many years there has been tension between modelers of stellar structure and 
cosmologists in that the modelers tended to derive ages for the Galaxy's oldest 
clusters that approached or exceeded the apparent age of the Universe. For much 
of that time the stellar astrophysicists had the advantage of better-constrained 
physics, but the recent resu lts from the WMAP mission lead to an age for the 
Universe of 13.7 ± 0.2 Gyr (iBennett et al. 20031 ). In other words, cosmologists 



can now claim to determine the age of the Universe to within about 1%. By 
contrast, the ages of individual nearby stars can be estimated to no better than 
10%, and that uncertainty is probably optimistic and does not fully take into 
account systematic effects. Stellar astrophysicists should see the WMAP results 
as a friendly challenge that we can meet and exceed, but only after significant 
exertion. We should be able to understand the closest stars at least as well as we 
can the entire cosmos! 



15 Sidebar 3: The Sun as a Benchmark 

There is exactly one stellar age that is both precise and accurate, that of the 
Sun, and it illustrates some of the inherent prob lems in determining ages. The 
bun IS 4, 567 ± 1 ± 5 Myr old (iGhaussidon 2007l V The extraordinary precision 



of 1 Myr represents measure ment error (individual measurements are precise to 
0.6 Myr ( Amelin et al. 20021 )). and the only slightly larger systematic error of 5 



Myr is due to uncertainty over the precise sequence of events in the early years of 
the solar system's history. That systematic error should lessen as we understand 
those events better. This age is determined from the decay of radionuclides. 

The Sun itself does not reveal its age in any of its observable properties. It 
is only because we can measure solar system material in the laboratory that we 
can establish the Sun's age with complete confidence; that is not possible for 
any other star. The Sun thus forms the only age benchmark available in all 
of stellar astrophysics, with all others being model-dependent to some extent. 
The extraordinarily precise age of the Universe from WMAP (13.7 ± 0.2 Gyr; 
Bennett et al. (20031 )) is based on some assumptions, but the accuracy is likely 



better than 10% at the very least. 

The Sun matters in another crucial respect, and that is for its absolute abun- 
dance scale. Abundances in stars are measured relative to the Sun, but the 
abundances that determine the structure of stellar models are absolute values. 
The recent controversy over CNO abundances in the Sun matters for ages because 
changing solar absolute abundances requires changes in all the models. 

For more on the Sun as a fundamental calibrator, see Christensen-Dalsgaard (20091 ) .[ 
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Figure 1: Normalized age probability distribution functions (PDFs) for isochrone 
placement of individual main sequence solar-type stars, chosen to show repre- 
sentative cases (adapted from Takeda et al. (20071 ) ). Stars C (brown) and D 
(red) show "well-defined" ages, having a single or at least a dominant peak. For 
star A (blue), an upper bound to age (2.9 Gyr) results plus a best-estimate age 
(1.2 Gyr). Only the lower-bound for star E (green) can be specified (12.2 Gyr). 
No meaningful age can be derived for star B (pink). Of the stars studied by 
iTakeda et al. f;2007l ). 2/3 had well-defined ages (in the sense they describe, with 
a well-delineated peak in the PDF) and of those only 25% had errors in age of 1 
Gyr or less. 
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Figure 2: A comparison of several sets of PMS models. For both 
panels the blue lines show the same mass tracks and isochrones from 
Siess. Dufour &: Forestini that are shown in Figure [3l except for 



a more limited range of m ass. Three sets of mode ls are shown: 
Siess. Dufour &: Forestini blue; iPalla &: Stabler (199^ ). green; and 



^ / I ■/ / (J 

D'Antona Sz Mazzitelli (19971 ). red. The three sets do not always cover the same 



masses or ages, and so a representative subset is shown. Mass tracks are shown 
on the left for four masses: 0.4, 0.6, 1.0, and 1.5 Mq. The cross near the bottom 
shows the effect of an uncertainty of ±100 K in T^g and ±0.1 dex in logLboi- 
Actual uncertainties in T^ff for PMS stars are larger (200-300 K) and also larger 
in Lboi (see Fig. [3]). The r ight panel shows isoch rones, marked with ages in Myr. 
The "0 Myr" isochrone of lPalla fc Stabler (19991 ) is their "birthline" where stars 
first become visible. Both sets of mass tracks and isochrones are broadly similar 
in many cases but deviate significantly at the youngest and oldest ages shown. 
The different nature of the mass tracks of iD'Antona &: Mazzitelli (19971 ) is es- 
pecially striking. Note the considerable uncertainty in age or mass that results 
from the models, independent of observational uncertainty. 



David R. Soderblom 



57 




Figure 3: An HR D for the Orion Nebula Cluster (ONC), adapted from 
Da Rio et al. f 20091 1. The ONC is one of the most-observed clusters of PMS stars 



and so illustrates the difficulties involved in studying the ages of very young low- 
mass stars (the highest-mass stars in the ON C are not shown) . Te mperatures were 
determined from spectral types (mostly from Hillenbrand (19971 )) to avoid the ef- 
fects on the photometry of reddening, veiling, accretion, and so on; this accounts 
for the quantization of most of the T^s values. Luminosities were determined 
from optical UBVRI photometry after correcting for reddening and veiling. The 
Rv value assumed for the reddening is 3.1, and the assumed distance to the ONC 
is 414 pc, as measured by radio interferometry. Individual values of reddening 
for each star were calculated using the multi-color photometry. The evol ution- 
ary tracks and isochrones shown are from Siess. Dufour &: Forestini (2000l ). The 
evolutionary tracks (dotted lines) have their masses shown along the 0.1 Myr 
isochrones, from 0.1 to 3.0 Mq. Several representative isochrones (solid lines) 
are plotted for 0.1, 0.5, 1, 2, 7, and 20 Myr, plus a ZAMS. The 1 and 7 Myr 
isochrones are bolder. The ages of these ONC stars scatter in Lboi enough to 
appear to range from ~ 1 — 10 M yr in age, peaking a t ~ 3 Myr. Despite this 
considerable spread in luminosity. Da Rio et al. (20091 ) did not conclude that a 
real age spread was clearly present because of the uncertainty added by stel- 
lar variability and scattered light from circumstellar material. Different sets of 
isochrones yield different results, bo th in absolute age and in the spread. For 
example, the evolutionary models of iPalla Stahler f 19991) lead to a mean a ge 
of about 2 Myr but with the age depending on mass (see Da Rio et al. (20091 )). 
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Figure 4: Observations of Prot in young clusters. The observed periods have 
been divided into the solar Carrington rotation (25.38 days) to yield an angular 
rate relative to the Sun. Shown are data for NGC 2264 (age approximately 5 
Myr), M50 (130 Myr); M34 (200 Myr); and the Hyades to show how rotation 
has dropped and converged by the Hyade s age of 625 Myr. Obse rvations of Prot 
are available for many other clusters (see Irwin &: Bouvier (20091 )) but these are 
representati ve of the trends. The da ta were provided by J. Irwin and the sources 
are listed in Irwin &: Bouvier f 20091 1. 
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Figure 5: Observations of the disk-integrated core of the Ca II K hne in the 
Sun. The observations are from a Lowell Observatory program that mo nitors 
the activity cycles of solar-type stars (see iHall. Lockwood &: Skiff (20071 ) for a 
description). The solar observations are fed into the same spectrograph used for 
the stellar work using an optical fiber. The time span of these data exceeds a 
solar cycle and the rise and fall of he overall activity of the Sun can be seen, 
together with variations. The scatter is entirely astrophysical, not noise. 
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Figu re 6: Observations of Ca ii H an d K activity at representative ages (adapted 
from iMamaiek fc Hillenbrand (200^ ) ) . The ordinate is a normalized index of HK 
strength, described in the text. Shown are observations for PMS groups (red; age 
5 Myr); the Pleiades (blue; 100 Myr); the Hyades (green; 625 Myr); and M67 
(black; 4 Gyr). The solar point is shown as a pink dot. Note in particular 
the range of HK indices seen in any one cluster and that the spread within the 
Pleiades, for example, exceeds the net change in the HK index in going from the 
mean Pleiades value to that of the Sun. 
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Figure 7: Observation s of lithium in clusters, using data from from 
ISestito Randich Shown are data for T Tauri stars at about 5 Myr 

(red), clusters from 30 to 50 Myr old (blue), the Pleiades at 100 Myr (green), 
and clusters at about 600 Myr (black) . The data are shown in observational units 
to illustrate the inherent scatter seen. This scatter is real and not error. Note 
that the T Tauri stars are undepleted in Li and exhibit little scatter in EW, while 
stars in the Pleiades have a spread in EW of 1 dex for {B — V)q> 1.0 (about KOV 
and later); the spread in abundance is even larger. 
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Figure 8: A comparison of ages determined from an HK activity in- 
dex and from isochrone placement. This sample consists of ~ 100 
solar- type star s that have detected plane tary companions. The HK 
ages are from Mamaiek fc Hillenbrand (20081 ) and the isochrone ages from 
Takeda et al. (2mi[ ). Note that the isochrone ages, on average, are about 1.5 



times the HK ages. 



